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Abstract 

We consider n non-intersecting Brownian motions with two fixed start- 
ing positions and two fixed ending positions in the large n fimit. We show 
that in case of 'large separation' between the endpoints, the particles are 
asymptotically distributed in two separate groups, with no interaction 
between them, as one would intuitively expect. We give a rigorous proof 
using the Riemann-Hilbert formalism. In the case of 'critical separation' 
between the endpoints we are led to a model Riemann-Hilbert problem 
associated to the Hastings-McLeod solution of the Painleve II equation. 
We show that the Painleve II equation also appears in the large n asymp- 
totics of the recurrence coefficients of the multiple Hermite polynomials 
that are associated with the Riemann-Hilbert problem. 

Keywords: non-intersecting Brownian motions, Riemann-Hilbert prob- 
lem, Deift-Zhou steepest descent analysis, Painleve II equation, multiple 
Hermite polynomials, recurrence coefficients. 

1 Introduction and statement of results 
1.1 Non-intersecting Brownian motion 

This paper deals with non-intersecting Brownian motions with two fixed starting 
positions ai, a2 and two fixed ending positions 61, &2- Let us first describe the 
general framework in which this paper fits. 

Let p,q G N. Consider sequences of real numbers {ak}^^i, {bi}1^i and 
sequences of positive integers {nk}^^i, {"t-;}^^]^ satisfying 

• Oi > • • • > Qp, 

• bi > ■ ■ ■ > bq, 
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Figure 1: Non- intersecting Brownian motions in case of one starting and ending 
point. Here the horizontal axis denotes the time, t G [0, 1], and for each fixed t 
the positions of the n non-intersecting Brownian motions at time t are denoted 
on the vertical line through t. 

Consider n one-dimensional Brownian motions (actually Brownian bridges) such 
that Uk of them start from the point at time t = 0, for k = 1, . . . ,p, and mi 
of them arrive in the point 5/ at time t — 1, for I — I, . . . , q, and such that the 
n particles are conditioned not to collide with each other in the time interval 

(0,1). ^ 

We are interested in the limiting behavior where the number of Brownian 
particles n tends to infinity in such a way that each of the fractions rik/n, k — 
1, . . . ,p and mi/n, I = 1, . . . , q has a limit in (0, 1). To obtain non-trivial limiting 
behavior we assume that the transition probability density of the Brownian 
motions scales as 

PN{t,.,y)^^e^^^^-y^' (1.1) 

V ZTTC 

where A'' is a parameter (the inverse of the overall variance) that increases as n 
increases. In a non-critical regime we will simply take N = n. 

Under these assumptions, it is expected that the particles will asymptotically 
fill a bounded region in the time-space plane (tx-plane). It is also expected 
that for every t G (0, 1), the particles at time t are asymptotically distributed 
according to some well-defined limiting distribution. See Figures [T] and [2] for 
possible behaviors when p — q — 1 and p — q = 2, respectively. 

It is natural to ask for an explicit description of this limiting distribution. 
Such a result is known in the classical case where p = q = 1 , which is connected 
to Dyson's Brownian motion [19j . Here it is known that up to a suitable scaling 
and translation, the Brownian particles at time t e (0, 1) have the same joint 
distribution as the eigenvalues of the Gaussian unitary ensemble (GUE) of size 
n X n, and the limiting distribution as n — > cx) is given by Wigner's semicircle 
law. More precisely, if all Brownian motions start from a := ai at time t = 



2 



and end up at b := bi at time t — 1, and if A'^ = n/T, then the non-mtersecting 
Brownian particles at time t G (0, 1) are asymptotically distributed on the 
interval [a(t), P(t)] with endpoints 



a(t) = (1 - t)a + tb- y/4Tt{l - t), (1.2) 
(3{t) = (1 - t)a + tb+ ^ATt{l - t), (1.3) 

and with limiting density of the particles given by the semicircle law on that 
interval 

1 



2TrTt{l - t) 



y/m)~x){x-ait)), x€[a{t),P{t)]. (1.4) 



By varying t G [0, 1], the endpoints (II. 2p . (II. 3|) parameterize an ellipse in the 
time-space plane, see Figured] 

Explicit descriptions of the limiting distribution of the non-intersecting Brow- 
nian motions are also known when p = 1, q > I. Also in this case there is an 
underlying random matrix model H El [23 EU US] ■ 

The limiting distribution is also known when p = q = 2, uii = rii = m2 = 
n2 = n/2, N = n, and (ai— a2)(6i — 62) < 2. In this case the limiting distribution 
as n — ^- CX3 of the Brownian particles at time t G (0, 1) is obtained from a certain 
algebraic curve of degree four 15j. In contrast to the previous cases, however, 
it is not known if the Brownian particles for finite n can be described as the 



eigenvalues of a random matrix ensemble. See Figure 2(b) for an illustration of 
this case. 

1.2 Separation of the endpoints 

The main results of this paper will concern non-intersecting Brownian motions 
with two starting points ai > 02 at time t ~ 0, two ending points 61 > 62 at 
time i = 1, and in addition 

ni = nil, n2 = m2, (for all n) (1-5) 

and if we put 

Pi = — , P2 = — (1-6) 

n n 

which are varying with n, then we assume that 

— = Pi = pI + 0{l/n), n^oo, (1.7) 



n 
n2 

n 



= P2 = p*2 + 0{l/n), n^oo, (1.8) 



for certain Hmiting values pl,P2 G (0, 1). We also assume that N increases with 
n such that 

T^n/N>0 (1.9) 

is fixed. 

If the assumptions (|1.5p - (|1.9p hold, and if the separation between the start- 
ing points oi and 02, and the ending points bi and &2 is large enough, then in 
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analogy with (|1.2p - (|1.4p . one would expect the Brownian particles to be asymp- 
totically distributed on two disjoint ellipses in the ia;-plane, whose intersections 
with the vertical line through t are given by the two intervals [a* , PI] and , 
with 



a* = a* (i) = (1 - t)aj + tbj - ^ Ap*Tt{l - t), 
/3* = P*{t) = (1 - t)a, + th, + ^4p*T<(l - t), 



(1.10) 
(1.11) 



for j — 1,2, and with limiting densities on these two intervals given by the 
semicircle laws 



^^y^(/3;-x)(x-a*), X £[«;,/?;], j = i,2. 



27rrt(l - 

for each t e (0, 1). This situation is illustrated in Figure 2(a) Note that 

r/9,' 



(1.12) 



1 



2^rt(i - 1) J,, 



(/3* - x){x - a*)dx = p*. 



We derive the precise condition for this two-ellipse scenario to happen. It is 
clear that a necessary condition is the disjointness of the two ellipses. 

Lemma 1.1. (Disjointness of the two ellipses) The two ellipses parameterized 
by p.lOp - p. lip are disjoint if and only if 



(ai - a2)(&i -b2)>T + v^)^ 



[1.13) 



PROOF. The two ellipses are disjoint if and only if ai{t) > /3|(t) for all 
t e (0, 1). From (frTOl) - (frTT|) this leads to the condition 



(1 - t){ai - 02) + t{bi - &2) > v/4rt(i - t) + , t e (0, 1), 



(1.14) 



which after putting u = \J^t is equivalent to 

1*2(61 - 62) - 2it\/r + ^) + (ai - 02) > 0, 



It £ (0, 00). 



The left-hand side is a quadratic expression in it, whose discriminant is negative 
if and only if (jl.l3p holds. The lemma then easily follows since 61 > 62 and 
ai > 02. □ 



We will call (|1.13p the case of large separation. 

Definition 1.2. (Large, critical and small separation) For each n, consider n 
non-intersecting Brownian motions with two starting points ai > a2 at time 
t = and two ending points 61 > 62 at time t — 1, and assume that the 
hypotheses (|1.5p - (|1.9p hold. If (jl.l3p holds, we say that we are in a situation 
of large separation of the endpoints. If instead 

{a,^a2)ibi-b2) <Ti^, + ^f, (1.15) 
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Figure 2: Non- intersecting Brownian motions with two starting and two ending 
positions in case of (a) large, (b) small, and (c) critical separation between the 
endpoints. Here the horizontal axis denotes the time, t G [0,1], and for each 
fixed t the positions of the n non-intersecting Brownian motions at time t are 
denoted on the vertical line through t. Note that for n — > oo the positions of 
the Brownian motions fill a prescribed region in the time-space plane, which is 
bounded by the boldface lines in the figures. Here we have chosen iV = n = 20 
and pi = P2 = 1/2 in each of the figures, and (a) ai = —02 = 1, bi ~ —62 = 0.7, 
(b) oi = —02 = 0.4, bi ~ —62 = 0.3. and (c) rii = —02 = 1, fei = —62 = 1/2, in 
the cases of large, small and critical separation, respectively. 
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we are in a situation of small separation, and if 



(ai - a2)(&i - 62) = Ti^pl + y^p*y, (1.16) 

we are in a situation of critical separation between the endpoints. In the latter 
case, we define the critical time <crit G (0, 1) as 

(fli - 02) + (Ol - O2) 

The non-intersecting Brownian motions corresponding to each of these three 



cases are illustrated in Figure 2(a) 2(c) In the case of critical separation, the 



time icrit is precisely the time where the two ellipses with endpoints parameter- 



ized by (|1.10p - (|l.lip are tangent, cf. Figure 2(c) 

We may equivalently view the critical behavior in terms of the temperature 
parameter T — n/N. With aj, bj, and p* fixed for j = 1,2, there is a critical 
temperature 

_ (ai - a2)(6i - 62) 

crit — 



such that T < Tdit (low temperature), T > Tcrit (high temperature), and T = 
Tcrit (critical temperature) correspond to large, small, and critical separation, 
respectively. 

1.3 Large separation: decoupling of the Brownian motions 

The case of small separation of the endpoints was considered in '15j for the 
special case when pi = P2 = 1/2. In the present paper, we will focus instead on 
the cases of large and critical separation. As is to be expected, in these regimes 
the Brownian motions asymptotically decouple into two separate groups, as is 



evidenced by Figures 2(a) and 2(c) This is our first main theorem. 



Theorem 1.3. (Decoupling of Brownian motions) Consider n non-intersecting 
Brownian motions with two starting points ai > 02 at time t — and two ending 
points bi > 62 ttt time t ^ I. Assume that the hypotheses (jl.5p - (ll.9p hold, and 
assume that either 

• there is large separation of the endpoints (jl.l3p . or 

• there is critical separation of the endpoints (|1.16p and t 7^ tcrit, where tcrit 
is given by (|1.17p . 

Then as n —> 00 the Brownian particles at time t G (0, 1) are asymptotically 
supported on the two disjoint intervals [a1{t) , I3'[{t)\ and [02(^)7 /52(^)] given by 
(|1.10p - (ll.lip . with limiting densities given by the semicircle laws (jl.l2p . 

In the case of critical separation, we strongly expect that the conclusion of 
the theorem should also be valid when t = icrit- However, we will not consider 
the critical time in this paper. 

Our proof of Theorem 11.31 follows from a steepest descent analysis of the 
underlying Ricmann-Hilbcrt problem that will be described in Section [1.41 The 
Riemann-Hilbert problem is of size 4x4 and it was also used in |I5j to analyze 
the case of small separation. In the (apparently) simpler case of large separation 
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one expects that for large n, the 4x4 RH problem asymptotically decouples 
into two smaller RH problems of size 2x2. This is indeed the case, but in 
order to show this, we need a preliminary transformation where we introduce 
auxiliary curves in the complex plane and subsequently perform a Gaussian 
elimination step in the jump matrix of the RH problem, serving to annihilate 
some undesired entries of this matrix. This Gaussian elimination step is similar 
to the so-called global opening of the lens discussed in [31 [H [TSl HHj ■ 

It follows from our analysis that the interaction between the two groups of 
Brownian particles decays exponentially with n in case of large separation, and 
polynomially (like a power n^^^^) in case of critical separation at a non-critical 
time. In the limit when n oo the particles will then indeed be distributed 
inside two disjoint ellipses in the time-space plane. 

1.4 Riemann-Hilbert problem 

The non-intersecting Brownian motions described in the previous subsections 
are related to the following Riemann-Hilbert problem (RH problem) which we 
already alluded to above. The RH problem was introduced in [14j as a general- 
ization of the RH problem for orthogonal polynomials in [22], see also [37]. In 
accordance with Section 11.11 we will state the RH problem for general numbers 
of starting and ending positions p, q of the Brownian motions, although in our 
applications we will eventually take p = q = 2. 
Define weight functions 



The RH problem consists in finding a matrix-valued function Y{z) = 
yni,...,«p;mi,...,m,(^) of sizc p + q hy p + q such that 

(1) Y{z) is analytic in C \ M; 

(2) For a; e M, it holds that 



where Ik denotes the identity matrix of size k; where W{x) denotes the 
rank-one matrix (outer product of two vectors) 



and where the notation Y-^{x),Y^{x) denotes the limit of Y{z) with z 
approaching x G M from the upper or lower half plane in C, respectively; 

(3) As z — + CX3, we have that 



w'i,fc(:E) = e-^(^'-2-'=-), fc = l,...,p, 
W2Ax)^e-^^'''-''"^\ 1 = 1,..., 



(1.18) 
(1.19) 




(1.20) 




(1.21) 



Y{z) = (/p+, + 0(l/z))diag(z' 




,...,z-™''). (1.22) 
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The RH problem has a unique solution |14j that can be described in terms 
of certain multiple orthogonal polynomials (actually multiple Hermite polyno- 
mials); details will be given in Section \TM 

Let us explain the connection between the non-intersecting Brownian mo- 
tions and the RH problem in the case p = q = 2. It is well-known [25] . 
see also [23J [5B , that the distribution of the non- intersecting Brownian mo- 
tions at time t e (0,1) describes a determinantal point process, determined 
by an associated correlation kernel. According to [H] the correlation kernel 
K{x,y) = Kni. ,712 'ifni ■,m2 

{x,y) can be expressed in terms of the solution to the 

RH problem as 



/ wi,i(a;)\ 


V ; 



(1.23) 

By general properties of determinantal point processes, Theorem 11.31 then 
comes down to the statement that under the conditions of Theorem ll.3[ the 
limit of —Kni.n2;ni,n2{Xi x) ss n ^ oo cxists and is equal to 



(1.24) 

for j = 1,2. This is what we will show in Sections [5] and [31 

Our method will also allow us to obtain the local scaling limits of the cor- 
relation kernel that are common in random matrix theory and related areas, 
namely the sine kernel in the bulk and the Airy kernel at the endpoints of the 
intervals [a^,/3j['] and [q!2,/?2]- ^^^^ ^'^^ go into details about this in this 
paper. 

1.5 Critical separation and the double scaling limit 

Assume that the endpoints are such that 

(ai-a2)(6i-62) = (v^+ V^)'. (1-25) 

We put 

N = n/T 

where T can be interpreted as a temperature variable. For T = 1 (critical 
temperature) we have by (|1.16p . (|1.25p and Theorem ll.3l that in the large n limit, 
the particles fill out two ellipses, which are tangent to each other at the critical 



time tcrit, see Figure 2(c) By varying T around the critical value Tcrit = 1 we 
move from a case of disjoint ellipses (for T < 1) to a case of small separation 
(for T > 1), where the two-ellipses scenario is not valid anymore. Hence we 
see a phase transition in the case of critical separation, which is clearly seen at 
the critical time tcrit- At a non-critical time t ^ tcrit the phase transition is less 
obvious, but there is also a nontrivial transitional effect. Indeed, the distribution 
of particles at time t ^ tcrit in the case of small separation differs from the 
distribution of two semicircle laws of two disjoint intervals, which it is for large 
separation. Hence the endpoints of the intervals do not depend analytically on 
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the starting and ending points, which indicates the phase transition. It is a 
surprising outcome of our analysis that the phase transition for the case where 
t 7^ tcrit can be described by the Painleve II equation. 

In the case of critical separation we investigate the behavior of the Brownian 
particles in a double scaling limit where n,N ^ oo, and simultaneously T — )■ 
Tcrit = 1- More precisely, we consider the endpoints ai, 02, 61, 62 fixed so that 
(|1.25p holds for all n. The temperature T = Tn = n/N is varying with n as 
follows: 

r„ = l + Ln-2/3, (1.26) 

where L is an arbitrary real constant. 

We will show that in the double scaling regime described above, the steep- 
est descent analysis of the Riemann-Hilbert problem leads in a natural way to 
a model Riemann-Hilbert problem related to the Painleve II equation. More 
precisely, we will be led to the construction of a local parametrix that can be 
mapped onto the model RH problem [5U] satisfied by the ^'-functions (Lax pair) 
associated with the Hastings-McLeod solution of the Painleve II equation 

q"is)^sqis) + 2qHs). (1.27) 

The Hasting-McLeod solution is the special solution q{s) of p.27p which is 
real for real s and satisfies q{s) ~ Ai(s) as s ^ 00, where Ai denotes the usual 
Airy function. The precise form of the model RH problem will be described in 
Section [3131 

The Hastings-McLeod solution of the Painleve II equation also appears in 
the famous Tracy- Widom distributions [331 ISS] for the largest eigenvalues of 
large random matrices. It also appears in the critical unitarily invariant matrix 
models, where the parameters in the model are fine-tuned so that the limiting 
mean eigenvalue density vanishes quadratically at an interior point of its support 
[6l[Tl] . In this case it leads to a new family of local scaling limits of the eigenvalue 
correlation kernel that involve the ^'-functions associated with q{s). 

In our situation the Painleve II equation does not manifest itself in the local 
scaling limits of the correlation kernel. The construction of the local parametrix 
is done at a point xq strictly outside of the support and it does not influence the 
local correlation functions for the positions of any of the particles. The point 
Xq does not seem to have any physical meaning. 

We emphasize that our asymptotic analysis will be only valid when t ^ tcrit ■ 
At the critical time t — <crit where the two ellipses are tangent, one is led 
to a considerably more difficult, multi-critical situation. Here one expects the 
appearance of a model RH problem related to some as yet unknown fourth 
order ODE. As already mentioned, we will not attempt to study this case in the 
present paper. 

1.6 Generalities on multiple orthogonal polynomials 

While the appearance of the Hastings-McLeod solution of the RH problem does 
not affect any of the local scaling limits, it is felt by the recurrence coefficients 
of the multiple Hermite polynomials. These polynomials appear in the solution 
of the RH problem given in Section 11.41 

To state the results, let us first recall some generalities on multiple orthogonal 
polynomials in the sense of [14,. In accordance with Sections 11.11 and 11.41 we 
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will again give the definitions for general values of p and g, although in our 
applications we will eventually take p = q = 2. 

Definition 1.4. (Multiple orthogonal polynomials; cf. [T3^) Let p,q G N he two 
positive integers. Let there be given 

• A (finite) sequence of positive integers ni , 712 , . . . , rip £ N; 

• A sequence of weight functions wi^i{x),wi,2{x), . . . ,wi^p{x) : K — > R; 

• A sequence of positive integers mi, TO2, . . . , G N; 

• A sequence of weight functions W2,i{x),W2,2{x), . . . ,W2,q{x) : R ^ M. 

Put n := {ni, . . . ,np), |n| := X]fc=i "-fc ^^^"^ similarly for m and |m|. Assume 
|n| = |m| + 1. We say that a sequence of polynomials Ai{x), ^2(2;), . . . , Ap{x) is 
multiple orthogonal with respect to the above data if (i) the polynomials Ak{x) 
have degrees bounded by Uk ~ 1: 

degAk<nk-l, k = l,...,p, (1.28) 

and (ii) the function 

p 

Q(a;) :=^Afc(x)u;i,fc(a;) (1.29) 
fc=i 

satisfies the orthogonality relations 

/oo 
Q{x)x^W2,i{x) dx^Q, (1.30) 

for j = 0, 1, . . . , m; — 1 and I — \, . . . ,q. 

Note that (|1.30p states that Q{x) has rai vanishing moments with respect 
to the weight W2,i{x), I — 1, . . . , q. 

A schematic illustration of Definition 11.41 is shown in Figure [3] Let us com- 
ment on this figure. The left part of the figure shows the polynomials Ak{x) 
and their corresponding number of free coefficients Uk, k — 1, . . . ,p. The mid- 
dle part of the figure shows how the polynomials Ak{x) should be assembled 
into the function Q{x) defined in p.29p . Finally, the right part of the figure 
schematically shows the orthogonality relations of Q{x) with respect to the dif- 
ferent weights W2^i{x), and it shows next to each weight W2.i{x) also the number 
of vanishing moments mi of Q{x) with respect to this weight. 

We will refer to the polynomials Ak{x) in Definition 11.41 as multiple or- 
thogonal polynomials (MOP). Note that these polynomials were called multiple 
orthogonal polynomials of mixed type in [T3] and mixed MOPs in 2 . We will 
also find it convenient to use the vectorial notation A(a:) — {Ai{x)^ . . . , Ap{x)). 
To stress the dependence on the multi-indices n, m we will sometimes write 
A(x) = An,m{x) and similarly Aj(x) — (^j)n,m(a;). 

The coefficients of the multiple orthogonal polynomials in Definition 11.41 
can be found from a homogeneous linear system with |n| unknowns (polynomial 
coefficients) and |m| equations (orthogonality conditions). The restriction |n| = 
|m| + 1 in Definition 11.41 guarantees that this system has a nontrivial solution. 
In fact, the solution space to this linear system will be at least one-dimensional. 
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ni Ai{x) 



-Wl.l 



W2^ 



mi 



n2 A2(x) ^Q(2.) W2,2 ^ ^2 



Up Ap{x) ^wi^y 



W2,q 



Figure 3: The figure shows a schematic ihustration of the muitipie orthogonal 
polynomials Ak{x) in Definition 11.41 



This corresponds to the fact that the MOP are only determined up to some 
multiplicative factor. If the MOP are unique up to a multiplicative factor then 
the pair of indices n, m is called normal [Tl] . 

The fact that the multiple orthogonal polynomials are only determined up 
to some multiplicative factor allows for different choices of normalization. 

Definition 1.5. (Normalization types of MOP; cf. [Mj) Assume the data in 
Definition 1 1 .41 and assume that n, m is a normal pair of indices. Then the MOP 
in Definition 11.41 are said to satisfy 

• the normalization of type I with respect to the Ith index, I € {1, . . . , g}, if 
the mjth moment of Q{x) with respect to W2.i{x) is equal to one, i.e., if 

Q{x)x""w2Ax) dx^l. (1.31) 



the normalization of type II with respect to the kth index, k G {1, . . . ,p}, 
if the leading coefficient of Ak{x) is equal to one, i.e., if 

Akix) = + 0(x"^--2). (1.32) 



The vectors of MOP An,rn(a;) corresponding to the above normalizations 
will be denoted as An^m (x) and A^^^i^^ (x) , respectively. 

The above normalizations might not always be possible. The type (/, I) 
normalization is not possible in those cases where the integral on the left side of 
p.3ip is equal to zero. Similarly, the type (//, k) normalization is not possible 
in those cases where the fcth polynomial Ak (x) has degree strictly smaller than 
Uk - 1. 

As mentioned above, the MOP appear in the solution of the RH problem in 
Section [LH Let us describe this is somewhat more detail. 

Recall that in Definition 1 1.41 we needed the condition |n| = |m| + 1 to ensure 
the existence of the multiple orthogonal polynomials. But let us now assume 
that |n| — |m|. In this case, the definition of MOP makes no sense. Indeed, the 
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coefRcients of the MOP would then solve a homogeneous linear system with as 
many equations as unknowns, which has in general only the trivial solution. 

Therefore, to apply the definition of MOP in a meaningful way for a pair 
of multi-indices satisfying |n| = |m|, we should first adapt the multi-indices. 
There are essentially p + q ways to proceed: 

1. One can increase one of the components n^, i.e., one can work with the 
pair of multi-indices n -I- e^, m, for some k g {1, 2, . . . ,p}. 

2. One can decrease one of the components m;, i.e., one can work with the 
pair of multi-indices n, m — e;, for some / G {1, 2, . . . , g}. 

Here denotes the vector which has all its entries equal to zero, except for the 
fcth entry which equals one. The length of should be clear from the context. 

Let us discuss the MOP corresponding to each of the p + q pairs of multi- 
indices above. In the case of a pair of multi-indices n-t-e^, m, we are dealing with 
multiple orthogonal polynomials where the fcth polynomial Ak, k E {1,2, . . . ,p} 
has increased degree; it will be natural to normalize the resulting MOP such 
that this fcth polynomial Ak is monic, i.e., to work with a normalization of type 
(//, fc). This leads to the vector of MOP 

AL+etU(^) = ((^l)n+etU(^)' ' ' ' ' i\)n+'Li^))- (1-33) 

In the case of a pair of multi-indices n, m — e; , we are dealing with multiple 
orthogonal polynomials where the Ith orthogonality condition, I e {1,2, ... ,q} 
has a decreased number of vanishing moments; it will be natural to normalize 
the resulting MOP such that this omitted moment equals one, i.e., to work with 
a normalization of type (/, I). This leads to the vector of MOP 

Ai';2-e,(^) - ((^i)L';^-e,(^),---,(^p)L';^-e,(^))- (1-34) 

Of course we are assuming here that all type I and type II normalizations 
in (jl.33p and p.34p exist. It turns out that the existence of each of these p + q 
vectors of MOP is equivalent to a single condition 15], to which we will loosely 
refer here as the solvability condition. This condition will always be satisfied in 
the case of Gaussian weight functions l|1.18p - (|1.19p . 

Now we can use the above MOP to solve the RH problem in Section 11.41 
Indeed, the p + q vectors of MOP in p.33p . (|1.34p are row vectors of length p, 
and they can therefore be stacked into the first p columns of a matrix of size 
{p + q) X {p + q). Denote with Y{z) such a matrix. The entries in the remaining 
q columns of Y{z) are defined as Cauchy transforms of the functions Q'^^^'^^ 
and Q^^'') defined as in (ll.29p . More precisely, the last q entries in the rows 
fc = 1, ... ,p of Y{z) are defined by 



2m J_^ X- z 

and in the rows p + k, . . . ,p + q by 

1 Q^'^''\X)W2,1{X) 



2ni 7-00 X — z 
The fact of the matter is the following. 



dx, I = 1, . . . ,q, (1.35) 



dx, I = 1, . . . ,q. (1.36) 
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Theorem 1.6. (Solution to the Riemann-Hilbert problem; cf. II4I) Let n, m 
with |n| = |m| be such that the solvability condition holds. Then there exists a 
unique solution Y{z) = l^n.mC-z) of the RH problem in Section \l.4\ This solution 
is given by ^ 

Y{z)^DY{z), 

where Y{z) is the matrix constructed from the MOP as described above, and 
where D — diag(/p, —2'KiIq). 

For example, in the special case where p = q = 2, the solution Y{z) of the 
Riemann-Hilbert problem in Section [1.41 is given by the 4x4 matrix 



fA \(II,2) \(II,2) 
Y(z) ^ D X \^^>n+e2,m l^2jn+e2,m * * 



(1.37) 



where D := diag(l, 1, — 27ri, — 27ri), and where the entries denoted with * are 
certain Cauchy transforms as in (fT35ll and (fO^ . The latter entries will be 
irrelevant in what follows. 

1.7 Recurrence relations for multiple Hermite polynomi- 
als 

In analogy with the three-term recurrence relations for classical orthogonal poly- 
nomials on the real line, one can show that the multiple orthogonal polynomials 
in Section [LH satisfy certain p + q + 1 term recurrence relations. We will state 
these relations in the case of multiple Hermite polynomials, i.e., when the weight 
functions of the MOP are given by the Gaussians ()1.18p - (|1.19|) . For simplicity, 
we assume throughout that p = q — 2, as in p.37p . 

Define the next term in the asymptotic expansion of Y{z) in p.22p as 

Y{z)= (^/+^ + (^^^^ diag(^"Sz"%z-™Sz-'"^). (1.38) 

The entries of the matrix Yi — (Yi)n.m in (|1.38p will be denoted by {cij)f j^^. 

Proposition 1.7. (Recurrence relations for multiple Hermite polynomials) As- 
sume p — q — 2 and let the weight functions be defined by (|1.18p ~ (|1.19p . Then 
the multiple Hermite polynomials in (|1.37p satisfy the 5-term recurrence rela- 
tions 

[11,1) 

31- 

Cl,2C2,l(^l)n+e2!m " Cl,3C3,l (^1 )n^m-ei " Cl,4C4,l (^1 )n^m-e2 ' (1-39) 



(^l)^+eV+ei,m+ei = [ ^ - {I - t)a, - tb, + (A, 

Cl,3 



(^l)n-Hei+ei,m+e2 = ( 2: - (1 - t)ai - tfea + ^ I (^l)n+ei,m 

- Ci,2C2,l(^l)i+e2m " Cl,3C3,l (^1 )n^m-ei " ^1,4044(^1)^^^^2^^, (1.40) 
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I A ^(-f-f.2) ( /, ,x I C2,lCi,3\ n(//,2) 

(^2)„+e2+e2,m+ei = ( ^ " (1 " " ^^1 + ^ I (-42)n+e2,m 

- C2,lCl,2(A2)i+e,^,m - C2,3C3,2(^2)n^m^2ei " C2,4C4,2 (^2 )n^m^2e2 ' (1-41) 

and 

(^2)„+e2+e2,m+e2 = U " (1 " i)a2 - ^^2 + ^ I (^2)n+e2,m 

- C2,lCi,2(A2)|/+e^^^_^ - C2,3C3,2(^2)n^m^2ei - C2,4C4,2 (^2 )n^m^2e2 • (1-42) 

Here we denote with Ci_j the entries of the matrix Yi = (1^1)11. m 'in (|1.38p . 

The proof of Proposition 1 1 . 71 will be given in a more general setting in Sec- 
tion 15.11 The explicit form of the first term in the right-hand side of each of 
(jl.39p - (|1.42p is only valid under the assumption of Gaussian weight functions 
(jl.l8p - (|1.19p : the explicit form of these terms will be established in Section [5.6.1l 

Note that the recurrence relations (|1.39p - p.42p contain several recurrence 
coefficients of the form CijCjj with i < j; it will be convenient to collect them 
in the 4 by 4 matrix 



fO C.i,2C2.l Ci,3C3,i Ci,4C4,i^ 

C2,3C3,2 C2,4C4,2 

C3,4C4.3 

\o y 



(1.43) 



It turns out that there exist certain connections between these recurrence coef- 
ficients. 

Proposition 1.8. (Relations between recurrence coefficients) Assume p ^ q ^ 
2 and let the weight functions be defined by (|1.18p - (|1.19p . Then the 2 by 2 
submatrix 

■Ci,3C3,i Ci,4C44 
,C2,3C3,2 C2,4C4,2y 



(1.44) 



of (jl.43p has row sums equal to t(l — t)j^ , k — 1,2, and column sums equal to 
t{l — t)^, 1 — 1,2. Next, assume that ni — mi and n2 — TO2- Then all the 
recurrence coefficients CijCj,i in (|1.43p can be expressed in terms o/ci_2C2,i and 
01^04^1 alone, by means of the following relations: 

C2,3C3,2 = Ci,4C4,l (1-45) 
Cl,3C3,l = *(1 - 0^ - Ci,4C4,l (1.46) 

C2,4C4.2 = <(1 - i)-^ - Ci,4C4,l (1-47) 
t^bi - b2fc3AC4,3 = (1 - tfiai - a2)'ci,2C2,l- (1.48) 

Note that p.45p - (|1.47p follow immediately from the stated row and column 
sum relations for the matrix (jl.44p . The latter will be established for general 
values of p and q in Section 15.41 see also Section 15.51 for a spectral curve in- 
terpretation of these relations. On the other hand, the equation (|1.48p will be 
established (in a slightly more general form) in Section [5.6.21 
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I. 8 Painleve II asymptotics for recurrence coefficients 

Finally we are in position to formulate the Painleve II asymptotics of the re- 
currence coefficients in p.39p - p.42p under the double scaling regime in Section 

II. 51 This is our second main result. We first consider the off-diagonal recurrence 
coefficients, i.e., the recurrence coefficients of the form CijCj^i with i < j. 

Theorem 1.9. (Asymptotics of off-diagonal recurrence coefficients) Assume 
the double scaling regime (|1.25p - (|1.26p . and let t e (0, 1) be a non-critical time, 
i.e., t ^ tcrit- Define the constants 

and 

s:=^-{plp*2)'^'{Vpl + Vpl)'^'LeR, (1.50) 
where L is defined in ()1.26|) . Then we have 

ci,2C2,i = -Kh^ibi - b2fq^is)n-^/^ + O(n-i), (1.51) 
Ci,4C4,i = Kh{l - t){ai - a2){bi - b2)q^ {s)n-^^^ + 0{n-^), (1.52) 

as n oo, where q{s) denotes the Hastings- McLeod solution to the Painleve 
II equation. The asymptotic behavior of the other recurrence coefficients CijCj^i 
with i < j is then determined by (|1.45p - (|1.48p in Proposition \1.8[ 



The key point of Theorem 1 1.9 1 is that it shows that the Painleve II equation 
shows up in the large n behavior of the recurrence coefficients in the case of 
critical separation at a non-critical time. 

Remark 1.10. (The case of large separation) Using the results in this paper, 
one can prove a similar result for the case of large separation of the endpoints 
(|1.13p . In this case, it can be shown that there exists a constant c > such that 

ci,2C2,i = 0(e-™), (1.53) 
ci,4C44 = 0(e-™), (1.54) 

as n — > oo. 

Remark 1.11. (The case of small separation) Performing a similar analysis of the 
results in [15., one can show that in case of small separation of the endpoints 
p.l5|) . and assuming the additional hypothesis pi = P2 = 1/2, the following 
expansions hold: 

ci,2C2,i = --77 7^(4 - (ai - a^fibi - 62)") + 0{n-^), (1.55) 

16(ai - 02)"^ 



ci,4C4,i - ^^V^(2 - (ai - «2)(6i - b2)) + O(n-i), (1.56) 



as n 



We have a similar theorem for the diagonal recurrence coefficients in (|1.39p - 
(|1.42p . i.e., for the first terms in the right-hand side of each of these equations. 
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Theorem 1.12. (Asymptotics of diagonal recurrence coefficients) Under the 
same assumptions as in Theorem \1.9\ we have that 



^ = -KVis)t(b, - 5.). M^l^i^^n-/^ + O(n-), (1.57) 

Cl,3 V Pi 



Cl,4 V fli — 02 



C2,3 V ai - 02 



= if^9^(.)t(6i - 52) J 1-^^) ^-2/3 + o(n-i), (1.60) 

C2,4 y P2 

as n ^ CX3. 



1.9 Phase diagram 

The main results of this paper and their relation to the other results known 
in the literature can be nicely summarized by means of a phase diagram. See 
Figure IH 

Let us comment on Figured Assume that the endpoints aj, bj, j ~ 1,2 are 
fixed and satisfy the critical separation (|1.25p . The horizontal axis in the figure 
denotes the time t G (0, 1) and the vertical axis denotes the temperature T > 0. 
The diagram is divided into different regions according to the behavior of the 
limiting distribution for n ^ oo of the non-intersecting Brownian motions at 
time t and temperature T. The region where T < 1 corresponds to the case of 
large separation; according to Theorem I 1 . 3l the limiting distribution is given here 
by two semicircles on the two disjoint intervals [a|(t), /3*(i)] and [a2(^)' /^l (Ol- 
At temperature T = 1, these two intervals meet each other at a certain critical 
time iciit- When T further increases, the intersection region between the two 
groups of Brownian motions starts to grow; this is indicated by the boldface 
curve in the middle of the picture. In the region below this curve but above 
T — I, the limiting distribution at time t is still supported on two disjoint 
intervals, but now with distribution described in terms of a certain algebraic 
curve of degree 4 [15], rather than semicircle laws. In the region above the 
curve, the limiting distribution is on one interval [15j . 

In case where p\ = P2 = 1/2; one can find an explicit description of the 
boundary curve in Figure U] from the results in [IS] ; it turns out that this curve 
is given by the equation 



T 

or equivalently 



oi — 02X^1 — < f bi — b2\^ t 



1-r 



y_ {ai-a2?{l-tf + {b,-b2?t^ _ ^^g^^ 



4t(l - t) 

The curve plotted in Figure [H corresponds to the choice of endpoints oi = —02 
61 = —62 = l/\/2 and hence tent = 1/2. 
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Figure 4: The figure shows the phase diagram for non-intersecting Brownian 
motions with two starting and ending points. We have chosen here p* = P2 = 
1/2, oi = —a2 = bi = —62 = l/-\/2, hence icrit = 1/2. The boldface curves 
denote the places where a phase transition takes place. The topmost curve has 
equation T = {t^ ^ t + l/2)/(<(l - i)), cf. (fmi) . 
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Figure 2] also displays the phase transitions between the different regions. 
On the horizontal line T — 1 the phase transition is described in terms of the 
Painleve II equation as shown in this paper. On the curve (|1.6ip one expects 
a description in terms of the Pearcey kernels; although for the case of two 
starting and two ending points this has not been strictly proven. For the case of 
non-intersecting Brownian motion with one starting and two ending points, the 
Pearcey kernels were obtained in [21 [ID] (in the equivalent setting of Gaussian 
random matrices with external source), see also [1] 13 [301 [35]. Finally, at the 
place where the two curves in Figure [U meet, one expects a phase transition in 
terms of an as yet unknown family of kernels. This is indicated by the question 
mark in the figure. 

1.10 Outline of the paper 

The remainder of this paper is organized as follows. In Sections [2] and [3] we 
apply the Deift-Zhou steepest descent analysis to the Riemann-Hilbert problem 
for multiple Hermite polynomials. We perform this analysis for the case of 
large separation in Section [2] and for the case of critical separation at a non- 
critical time in Section [3] leading to the proof of Theorem 11.31 In the critical 
case we are also led to a local parametrix for the RH problem in terms of 
the Painleve II equation. Next, we investigate the large n asymptotics of the 
recurrence coefficients of the multiple Hermite polynomials under the double 
scaling regime. This is done in Section [H leading to the proof of Theorems 11.91 
and 11.121 Finally, the Propositions 11.71 and 11.81 are established in Section [3 

2 Steepest descent analysis in the case of large 
separation 

In this section we analyse the non-intersecting Brownian motions in case of 
large separation between the endpoints (|1.13p . Using the Deift/Zhou steepest 
descent analysis of the Riemann-Hilbert problem, we show that the interaction 
between the two groups of Brownian particles decays exponentially with n, 
thereby establishing Theorem ll.3l 

2.1 Starting RH problem 

Our starting point is the RH problem (|1.20p - (|1.22p with p — q = 2 and in 
addition ni — nii and n2 = m2. As in (|1.6p we write pj — rij/n. Without loss 
of generality we take N — n (i.e., T — 1). For T — 1 the case of large separation 
corresponds to (ai — a2)(6i — 62) > {\/Pi + Vpf)^- Since pj ^ p* as n ^ 00, 
we already assume that n is so large that 



(ai - a2)(6i - 62) > (VPi + Vp^)^- 



Thus Y satisfies the following RH problem. 

(1) Y{z) is analytic in C \ M; 

(2) For a; e M we have 
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with the rank-one block W given by 

W ^ [ (W2,l W2,2) = • (2.1) 

\Wi^2j ^ ' V^l. 2^2,1 ^1,2^2,2/ 

(3) As z — > oo, we have that 

Y{z) = (/4 + 0(l/z))diag(z"Sz"%z-"Sz-"^). 

The entries of the rank-one block W = W{x) in (|2.ip can be written explicitly 
as 

wi^k[x)W2,i[x) ^ e 2t(i-t)^ (-2.2) 

for fc, Z e {1, 2}. Recall that = n. 

It will be convenient to write the diagonal entries of (|2.ip as 

wi,i{x)w2,iix) = e-"fi^i(^) = e-"i^i(^\ (2.3) 
ii'i,2(a;)w'2,2(a;) = e-"f^^=^(") - (2.4) 

with pj — Uj/n, for j — 1, 2, and 

V^i(^) 2pi^(l - ^) - 2 ((1 - t)ai + i6i) a:), (2.5) 
^^^""^ '-^ 2p2t{l - t) ^""^ - 2 ((1 - t)a2 + tb2) x). (2.6) 

Our goal is to show that in the large n limit the 4x4 matrix valued RH 
problem for Y{z) essentially decouples into two smaller 2x2 problems with 
weight functions (12. 3p and (12. 4p . To show that this decoupling indeed occurs, 
we need to show that in some sense, the off-diagonal entries in (|2.ip can be 
neglected with respect to the diagonal entries of (|2.ip . More precisely, one 
expects that 

• around the interval (a*, (31), the (1,1) entry of (j2.ip is dominant (as n 
oo) with respect to the other three entries. 

• around the interval (a2,/3|), the (2,2) entry of (|2.ip is dominant with 
respect to the other three entries. 

Here a*, (3* for j = 1, 2 are as in (frTOl) - (frTT|) . 

Remarkably, these expectations are not confirmed by a straightforward steep- 
est descent analysis in which the first transformation of the RH problem is based 
on the two semicircle densities (jl.l2l) and the corresponding g-functions. This 
approach turns out to be successful only for t near the critical time ^crit defined 
in (|1.17p . When t is sufficiently close to however, one runs into difficulties 
since then 

• the (1,2) entry of (|2.ip blows up (i.e., becomes exponentially large when 
n — > oo) somewhere in the interval (a*, PI), and 

• the (2, 1) entry of (|2.ip blows up somewhere in the interval (a2,/32)- 
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Similar problems occur when t is close to 1, but then with the roles of the (1, 2) 
and (2,1) entries of W reversed. 

In order to prevent the blow-up of undesired entries, we make a first prelim- 
inary transformation that is described in the next subsection. The transforma- 
tion is different for the two cases < t < icrit and tcrit < t < 1. For definiteness 
we assume from now on 

o<t< i,Ht = . — ,y (2.7) 

(ai - 02) -I- [bi - 02) 

The case where icrit < t < 1 is similar and the corresponding modifications will 
be briefly commented on later. 

2.2 First transformation: Gaussian elimination in the jump 
matrix 

The first transformation of the RH problem is a Gaussian elimination step for 
the jump matrix, serving to annihilate some of the undesired entries in (|2.ip . 
This elimination step will be at the price of introducing new jump matrices on 
certain contours Fi, F2 in the complex plane. This transformation is similar to 
the so-called global opening of the lens discussed in [3] 0] [HI |^ . 
Similar to ([LTOl) - (frTT|) we define 

= aj{t) = (1 - t)aj + thj - ^ 'ipjt{\ - t), (2.8) 
= = (1 - Oflj + + \J^P]t{^ - t), (2-9) 

for j = 1,2. Since pj is varying with n, the quantities (|2.8p ^ ()2.9p are also varying 
with n. For n ^ 00 they tend to a* and (i* given by (|1.10p - (|l.lip with T = 1. 
Define the corresponding semicircle laws 

^{fij-x){x-aj), xe[a,,/3,], j = l,2, (2.10) 



2TTt{l - t) 



which are also (slightly) varying with n. 

We take a reference point a;o G (/32, cti) and we choose unbounded contours 
Fi, F2 in the complex plane, crossing the real axis in points xi, X2 so that 

P2 < X2 < xo < xi < ai. 

The contour Fi stretches out to infinity in the right half-plane (i.e., Rez ^ 00 
as 2 ^ 00 on Fi), and F2 stretches out to infinity in the left half-plane. The 
precise way to choose xq and the curves Fi, F2 will be described later. We orient 
these curves in the upward direction as in Figure [51 We then define a new 4x4 
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Figure 5: The figure shows the curves Fi, F2, the intervals (—00,2:2), {x2,Xi), 
(cci, 00), and the jump matrices on each of these five curves in the RH problem 
for X = X{z). 



matrix- valued function X — X{z) by 



X = Y 



X = Y 



X = Y, 
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W2 
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t V 



to the right of Fi, 



to the left of F2, 



between Fi and F2. 



(2.11) 

(2.12) 
(2.13) 



The matrix function X satisfies a new RH problem, with jumps on the contour 
RUFi UF2. The jump matrices are different on each of the five pieces (—00, X2), 
{x2tXi), (xi,oo), Fi and F2. They are shown in Figure O 
Thus X satisfies the following RH problem. 

(1) X{z) is analytic in C \ (R U Fi U F2); 

(2) On R U Fi U F2 we have that X+ = X^Jx with jump matrices Jx as 
shown in Figure O 



(3) As z ^ 00, we have that 

X{z) = (/4 + 0(l/z))diag(z"\z"%z 



-"1 2^"2^ 



(2.14) 
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For the asymptotic condition (|2.14p we note that 

' (Oi — 62)2 



exp 



1 - t 



(2.15) 



^2,1(2) 

Since 61 > 62 we have that (|2.15p tends to exponentially fast as Rez — > cxd, 
and so in particular as z ^ 00 to the right of Fi. Similarly, the inverse of ()2.15|) 
tends to exponentially fast to the left of F2. Thus the transformation (j2.11|) . 
(|2.12p . (|2.13p leads to the same asymptotic behavior for X as we had for Y . 

Note also that the jump matrices on Fi and F2 tend to the identity matrix 
as z — > 00. 

What we have gained is that in the jump matrices on the intervals (—00, 2:2) 
and (xi, 00), the first and second columns of the top right 2x2 block (|2.ip of 
the jump matrix have been eliminated, respectively. 

Remark 2.1. The description of the Gaussian elimination step above has been 
done under the assumption (|2.7p . i.e., < t < fcrit- The case where tcrit < t < 1 
is different since it requires other entries of the jump matrix to be eliminated in 
the different regions of the complex plane. In order to do so one would define 
X differently as (compare with ([^T^ ) 



X = Y 



X ^Y 



X = Y. 



( 1 












^ 1 








Wl. 


1 












1 





V 








1/ 


/I 







'\ 





1 














1 





Vo 








1/ 



to the right of Fi 



to the left of F, 



in the region between Fi and F2. 



The further steps in the steepest descent analysis will then be similar to the 
case where Q < t < icrit and we will not discuss this any further. 

Alternatively, the results for tcrit < t < 1 could be reduced to those for 
< i < icrit by virtue of the involution symmetry to be discussed in Section 
15.3.3) a further note about this will be given in Remark [ 



2.3 Second transformation: (^-functions 

The second transformation is the normalization of the RH problem at infinity. 
To this end we use the so-called g-functions. 

As said before, in the limit n — > 00 we expect the Brownian particles to be 
distributed on two separate intervals [aj;, /9*] and [q!2, /32] (recall t is fixed), with 
limiting densities given by the two Wigner semicircle laws (|1.12p . 

For finite n, we have defined aj and (ij by p.8p - p.9p and the densities ()2.10p 
that are varying with n. We use the g-functions corresponding to these varying 
semicircle densities. We define 

log(z ~ x)\/ (/3i — x){x ~ ai) dx, (2-16) 

/Ql 

\og{z - x)^/{P2 -x){x-a2) dx, (2.17) 

f Oil 



92{z) 



1 

2Tlt(l - t) 

1 

2TTt{\ - t) 
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where the logarithms are chosen with a branch cut on the positive and on the 
negative real axis, respectively. Hence gi (z) is an analytic function on C\ [ai , oo) 
while g2{z) is analytic on C \ (—00, /32]- Note that 



9jiz) =Pjlogz + 0(l/z), 
It follows by contour integration that 



00. 



1 



2t{l - t) 
1 

2t{l - t) 



ai + I3i 
2 

012 + P2 



((z-ai)(z-/3i))i/2 
{{z ~ a2){z ~ (32))'/' 



(2.18) 
(2.19) 



for all z in the domain of analyticity of gi and g2 ■ Here the branches of the square 
roots arc taken which are defined in C \ [q!i,/3i] and C \ [a2,(32], respectively, 
and which behave as z when z — > 00. 

On the interval [ai, /3i], the square root in ()2.18p has two distinct boundary 
values, one from the upper and one from the lower half plane of C. These bound- 
ary values are complex conjugate and purely imaginary. Similar statements hold 
for the behavior of the square root in (|2.19p on the interval [a2 , (32] ■ 

It follows from the above observations and from the definitions (|2.5p . (|2.6p . 
(frTUll - (frTT|) that 



9i. 



-{x)+9[,-ix) 



1 



92,+ i^) + 92,-i^) 



ta-t) 

PiVlix), 
1 

P2V^{x), 



ai + (3i 



a2 + (32 



xe[ai,/?i], (2.20) 



xe[a2,(32]. (2.21) 

After integration, we see that there exist constants li, I2 such that 

9i,+ + 9u- ~ PiVi ^ h on[ai,/3i], (2.22) 

52,+ + 52,- - = ^2 on[Q!2,/32]- (2.23) 



Moreover, from the signs of the square roots in (j2.18p and (|2.19p we have 



91, + + 91,- - PiVi < h onR\[ai,/3i] 

92, + + g2- - P2V2 < h onR\[a2,/32] 



(2.24) 
(2.25) 



The above equations (|2.22p - (|2.25p are the Euler-Lagrange variational con- 
ditions for the equilibrium measure under an external field. Indeed, the Wigner 
semicircle laws ()2.10p (after normalization by a factor l/pj) are the equilibrium 
measures in the presence of the quadratic external fields Vi, V2, respectively, see 
[roil33j . The fact that Vi in (|2.5p has a factor pi in its denominator can be inter- 
preted by noting that for pi — > 0, the external field Vi gets stronger and stronger 
and hence the (potential theoretic) electrostatic particles are pushed together 
onto a narrower and narrower interval, which is confirmed by ()2.8p - (|2.9p . 

Now we use the (^-functions to normalize the RH problem at infinity. We 
define a new 4x4 matrix-valued function T = T{z) by 



(2.26) 
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where 



G = diag (e-fi , e"^^ , e^' , e^^) , (2.27) 
L = diag g-i2/2+.^ _ (2.28) 

Here gi, g2 are the g-functions, Zi, ^2 are the variational constants defined in 
(|2.22p and (|2.23p . and k e C is a constant to be determined. 

The RH problem for T = T{z) is normalized at infinity in the sense that 

T{z) = h + 0{\/z) (2.29) 

as z — > oo. This follows from the fact that gj{z) — pj logz + 0(l/z) as z — > oo 
for j = 1, 2, so that 

^ng,{z) ^ ^np, (^ _^ 0{l/z)) = z"^ (l + 0(l/z)) 

as z ^ OO. Hence the factor G" in (|2.26p cancels the powers of z appearing in 
()2.14p . Also note that the similarity relation with the matrix L in (|2.26p does 
not change the asymptotics (|2.29p . 

The jumps for T{z) are on the same contours M U Fi U r2 as those for 
X{z), but with different jump matrices. To state the resulting jumps, we find 
it convenient to work with the following functions ('A-functions') 

Ai(z) :=-gi(z) + ^ + i(z2-2aiz), 

H^) :=-52(^) + | + «+ij(z'-2a2z), 

A3(z) ■.^g,(z)-L-^(z^-2hz), 

>^4{z) := g2(^:) - Y + K - 2(i^(-z^ - 2fe2z). 

Using (|2.30p and after a little calculation, we find that we can write all jump 
matrices in the RH problem for T{z) as shown in Figure [6l On {x2,xi) the 
jump matrix is 

(o^ I) with W ^ (^^iTx!^) {e-'^^'^ e"^^W). (2.31) 

Thus T satisfies the following RH problem. 

(1) T(z) is analytic in C \ (R U Fi U F2); 

(2) On R U Fi U F2 we have that r_|_ — r_ Jt with jump matrices Jt as shown 
in Figure [S) 

(3) As z ^ 00, we have that 

r(z) = /4 + 0(l/z). (2.32) 

2.4 Asymptotic behavior of the jump matrices 

In this subsection, we investigate the asymptotic behavior in the large n limit 
of each of the jump matrices in the RH problem for T — T{z) in Figure [51 Our 
goal is to show that all the jump matrices are exponentially close (as n — > 00) 
to the identity matrix, except for the jump matrices on the intervals (ai,/3i) 
and (a2,/32) which have an oscillatory behavior. To be able do so we still have 
the freedom to take the reference point xq and the contours Fi and F2 in an 
appropriate way. We also have the constant k in (|2.28p at our disposal. 
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Figure 6: The figure shows the jump matrices in the RH problem for T = T{z). 
Here W is defined in (^31]) . and \i{z) , . . . , \i{ z) are defined in (|2.30p . 



2.4.1 Preliminaries on ^-functions 

In what follows wc will make extensive use of the ^-functions which are essen- 
tially the derivatives of the A-functions. We define 



C,(z)-A;.(z)-|^^z, j = 1,...,4. (2.33) 



From p.30p we then see that 



(2.34) 



Inserting (|2.18p - (|2.19p and (I2.8p ~ (|2.9p into these expressions, we obtain after 
some simplification, 



= 2t(h) (-(1 - *)«2 + + ((Z - a2)[z - /32))l/2) 
Uz) = 2t(T^ (-(1 - *)«2 + - {{Z - «2)(Z - 132)Y'^) 



(2.35) 



From either ()2.34p or (I2.35P we see that is defined and analytic in the cut 
plane C \ [a^, /3fc] (with fc = j or fc = j — 2), but Jie^j{z) is well-defined on the 
cut as well. From (|2.35p we easily deduce the following behavior on the real 
line, see also Figure [71 
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Figure 7: The figure shows the functions ^i(x), (.^{x) for x G R \ [q:i,/3i] and 
^2(2;), £,4{x) for a; e M \ [0:2, /32], together with the two horizontal hne segments 
Re(a(x)) = Re(e3(x)) = Si for x e [ai,/Ji] and Re(6(a;)) = Re (^4(2;)) = S2 
for X G [0:2. /32]- Wc have here pi = p2 = 1/2, t = 1/3, — a2 = ai = 1, 
—62 = bi = 1 and (hence) ai « 0.33, /3i « 1.66, a2 ~ —1.66, /32 ~ —0.33 and 
xo = 0. The graphs of ^2 and ^3 intersect at the a;- value — \/2/6 « —0.23 and 
those of ^1 and ^4 intersect at the a;-value y/2/6 w 0.23. 

Lemma 2.2. (a) We have that Re^i and Re ^3 are constant on the interval 

Rea(a;) = Re 6 (a;) = Si := ^^^^ (-(1 - t)ai +th), 

for X e [q!i,/?i]. On R \ [q!i,/3i] we have that ^1 and ^3 are real, ^1 is 
strictly increasing, ^3 «s strictly decreasing, and 

l{Ux)+Ux)) = ^u xeR\[aul3i]. 

(b) VFe /laue that Re ^2 Re ^4 are constant on the interval [q;2,/?2]- 

Re6(a^) = ReUix) = ^2 := ^^^^—^ ("(1 " 0«2 +^62) , 

for X e [a2,/32]- On R\ [a2,P2] we have that ^2 and ^4 are rea/, ^2 is 
strictly increasing, ^4 is strictly decreasing, and 

\ (6W + = S2, ar e M \ [a2,/32]. 
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(c) We have 

E2 > Si. 

PROOF. All statements immediately follow from (j2.35p . except maybe the 
statement of part (c). Part (c) follows from the fact that 

~ ~ ^ ■((l-i)(ai-a2)-i(6i-&2)) 



' 2t{l - t) 

which is indeed non-negative because of our assumption ()2.7p . □ 



The main property of the ^-functions is contained in the following lemma. 
It will only be valid under the large separation assumption, as we will see in the 
proof. 

Lemma 2.3. There is a value xq € (/32,cn) such that 

^xq) > 6(a^o) > Uxq) > ^i{xo). (2.36) 

PROOF. We prove the lemma in two steps. In the first step we show that 
there exists xq G (/32,ai) such that 

6(a;o) > Uixo)- 
This we can do by giving xq the explicit value 



'"'~Vp^+Vp~2 2 +vpr+vpJ 2 ^"-'^^ 

so that clearly 

<xo<^— 

which means that xq lies between the midpoints of the two intervals, and so in 
particular a2 < xq < Pi. 

We have by ^M-^M and (|07|) that 



{ai - Xo){Pl - Xq) 
Xq 



ai + Pi \^ fPi — cti^'^ 



2 

Pi f cti + Pi - a2 - P2^ 



4pii(l - t) 



{Vp'i + Vp^YK 2 

= , (((1 - t){ai - a2) + t{bi - b2)f - 4t(l - t){^i + . 

(2.38) 

Because we are assuming that we are in a situation of large separation we have 
{\/Pi + \/P2Y < ('^1 ~ o,2){bi — ^2) SO that we obtain from (|2.38p 

(ai - xii){Pi - X(i) 

> . (((1 - t){ai - 02) + t{bi - b2)f - 4i(l - t){ai - a2){bi - 62)) 

{{I^t){ai-a2)-t{bi-b2)f. (2.39) 



(%/pr+ VP2) 
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Similarly, 

(2:0 - 0:2) (xo - 132) 



(((1 - - 02) + t{h - b^ jf - 4t(l - t){^+^)^) 



It then follows that both factors on the left-hand sides of ()2.39|) and (|2.40p are 
positive, so that 

P2 <xo < ai, 

as claimed. 

Furthermore, we obtain from p.39p and (I2.40p 



a/ (ai - xo){pi - xo) + \J (a;o - a2)(a;o - 

> |(l-i)(ai-a2)-i(&i-&2)|. (2.41) 

Since by (|2.35p we have 
2t(l-t) (6(2^0) -^4(2:0)) 



= \J {a\ - a;o)(/3i - xq) + V (2^0 - a2)(a;o - /32) 

-((l-t)(ai-a2)-t(&i-62)), (2.42) 

it follows from (|2.4ip that ^3(2:0) > ^4(2^0), as claimed as well. 

In the second step we prove that the string of inequalities (|2.36p holds for 
some value xq g (/32,Q;i) which, however, could be different from x^ used in 
the first step. We distinguish three cases, depending on whether ^3(/32) < S2, 
iz[f32) = S2, or C3(/32) > S2 where S2 is the value of Re ^2 (a;) = Re ^4(0;) for 
X e [a2,/32], see part (b) of Lemma 1^?^ 

If ^3(/32) < S2, then ^3(/32) < ^^{fi-i)-, and since by what we already proved, 
the inequality ^3(2:) > ^4(0;) holds for some x G (/32,Q;i), there also exists 
X* S (/32,ai) where equality ^3(2:*) = ^^ix*) holds. Then by the monotonicity 
properties of the ^-functions (stated in parts (a) and (b) of Lemma 12. 2p 

i2{x*)>i^{x*)^i^{x*)>i^{x*). 

Then by slightly increasing x* we find xq such that the inequalities (|2.36p hold 
(with strict inequalities). 

If ^3(/32) — ^2, then as above we have 

6(/32)=C3(/32)=e4(/32)>6(/32). 

Then for xq slightly larger than /32, we have (j2.36p with strict inequalities, since 
^2(^0) — +CXD and ^4(0:0) = —00, see formulas (|2.35p . 

If iz^li-i) > ^2, then ^3(^2) > £,2{P2)- Since ^2 is increasing we find £,2{ai) > 
£.2(^2) = ^2 > Si = '^3(0^1)1 see Lemma Thus there exists xq & (/32,cn) 
with 6(2^0) ^sixo)- We have ^4(^0) < 6(2:0), 6(2^0) < 6(2^0), and 

6(2^0) +6(2^0) = 2Si < 2S2 =6(2^0) +6(2:0). 
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again by Lemma [^?^ so that ^i(a;o) < ^4(2:0) and the inequahties (|2.36p hold. □ 

It follows from the proof that in most cases we may assume that the in- 
equalities (j2.36p are strict. Only if t = tcrit then Si = S2 and then we can only 
obtain 

6(2^0) = 6(a;o) > £,4{xo) = Ci{xo). 
2.4.2 Jump matrices on Fi and r2 

After the preliminaries on the ^-functions we consider the jump matrices on Fi 
and We want that the jump matrices on these curves, shown in FigureEl are 
exponentially close to the identity matrix as n ^ 00. Thus we want to choose 
Fi and F2 so that 

for certain constants ci, C2 that do not depend on n. 

We take the point xq satisfying the inequalities (|2.36p of Lemma 12.31 We 
have by ((T^ 

Mz) - Hz) - 92{z) - I + « - + I - ^i-r^^. (2.44) 

The constant k is still at our disposal. We choose it here so that ()2.44p vanishes 
for z = Xq. Thus 

A4(a;o) - A3(a;o) = 0. 



Since the derivative is 

= U{xo) - 6(2^0) < 0, 



±iX4x)-X,{x)) 



see (|2.33p . we can then find xi and X2 sufficiently close to xq so that P2 < X2 < 
Xq < xi < ai and 

A4(a;i) - A3(xi) < < A4(x2) - A3(x2). 

We will choose Fi and F2 so that they cross the real axis in xi and X2, 
respectively. It remains to show that they can be extended to infinity, while 
remaining in the open sets 

f^i := {z e C I Re (A4(z) - A3(z)) < 0}, (2.45) 
:= {z e C I Re (X^iz) - X^iz)) > 0}, (2.46) 

respectively. Note that fli and fl2 are open subsets of C, symmetric with respect 
to the real axis and such that 

Xi G fii, X2 S il2- 

The fact that we can indeed choose Fi and F2 this way, follows from the following 
lemma. 
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Lemma 2.4. For j = 1,2, let Q° denote the connected component of flj that 
contains Xj . Then 51° is unbounded. In addition, for each e > there exists 
R > so that 

{zeC\\z\> R, -Tr/2 + e < argz < 7r/2 - e} C Q", (2.47) 



{z e C I |z| > R, 7r/2 + e < argz < 37r/2 - e} C n°. (2.48) 

PROOF. Wc will establish these properties by using the maximum princi- 
ple for subharmonic functions, where we use that Regj as given by (j2.16p is 
subharmonic on C and harmonic on C \ [aj, f3j] for j = 1, 2. 

Suppose that fij (the closure of ^1°) has nonempty intersection with the 
interval [a2,/32]- Since is connected, and symmetric with respect to the real 
axis, it will then surround the point X2, so that Q2 niust be bounded and 

n[ai,/?i] =0. 

Then Regi is harmonic on Re 172 is subharmonic and therefore by (12.44p 

Re(A4- A3)(z) ^ Re g2iz) - Re gi{z) - ^i—y Re z + const, (2.49) 

is subharmonic on fij- By definition we have that Re (A4 — A3) > on 
with equality on Qfij- This is a contradiction with the maximum principle for 
subharmonic functions and it follows that 

n?n[a2,/32] = 0. 

Then if Q° is bounded we obtain a contradiction with the maximum principle 
in the same way. Thus fij is unbounded, and likewise Q2 unbounded as well. 
From (|2.49p we further see that 

Re (A4(z) - A3(z)) = -^j^Re z + 0(log(|z|)) (2.50) 

as z ^ 00, since gj{z) = pj logz + 0(l/z) as z — )■ 00. In addition, we have 

^(A4(z) - A3(z)) = U^) - 6(^) = + Oil/z) (2.51) 

as z ^ 00. From (|2.50p - (|2.5ip it easily follows that for |Imz| large enough, we 
have that Re(A4(z) — A3(z)) monotonically decreases as Rez increases. This 
implies that the domains f^i and U2 given in (|2.45p - (|2.46p both have only one 
unbounded component, which then coincide with Q° and since these are 
unbounded. It now also follows from (f^SUll that (ITTtI) and hold. □ 



Re (A4 — A3 



We conclude from Lemma 12.41 that the curves Fi, F2 can be extended to 
infinity so that 

< -ci < 0, on Fi, 
> C2 > 0, on F2 

for some constants ci, C2 > 0. Since as n — > cxd we are in a non-critical situation, 
the contours F^ and the constants Cj can be taken independently of n for n large 
enough. Then the off-diagonal entries in the jump matrices on Fi and F2 are 
uniformly exponentially small as n — s- 00, as required. 
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2.4.3 Jump matrices on M 

We next investigate the jump matrices in the RH problem for T on the various 
parts of the real line. Recall that the jump matrices are given in Figure E] and 
so we are interested in the behavior of Xk.+ — ^j,- on the real line for various 
combinations of j and k. 
We prove: 

Lemma 2.5. The following hold. 

(a) For j — 1,2, we have 

Aj+2.+ - Aj,_ =0 on [aj,l3,], (2.52) 
Xj+2,+ - Aj,_ < onR \ [aj,f3j]. (2.53) 



(b) We have 



Re(A4,+ - Ai) < on {-oo,xi), (2-54) 
Re (A3,+ - A2) < on (2:2,00). (2.55) 



PROOF, (a) From the definitions ([2?30| and l[23 1) -([2?6 l) it follows that for 

Aj+2,+ - Aj,_ = - gj- - Ij ~ pjVj, 

so that part (a) is a restatement of the Euler-Lagrange conditions l|2.22p - (|2.25p . 

(b) Recall that the constant k was chosen so that A3(a::o) = A4(a;o). Then in 
view of part (a) we know that 

A4(a;o) - Ai(a;o) = A3(a::o) - Ai(a::o) < 0, 
A3(a;o) - A2(a;o) = A4(a::o) - A2(a::o) < 0. 

Taking xi and X2 sufhciently close to xq (which we can do without loss of 
generality), we then have the inequalities (|2.54p - (|2.55p on the interval {x2,xi). 

From Lemma 12.21 we know that Re ^4 is decreasing while Re ^1 is strictly 
increasing on (— 00, ai]. Then we have for x < xq, 

ReiU - > Rei^i - ^i){xo) = ^4(2:0) - U^o) > 

where the last inequality holds because of Lemma 12.31 Since 

^Re(A4-Ai) = Re(C4-6) 

it then follows that Re (A4 — Ai) is strictly increasing on (— oo,xo). Since we 
already proved that the inequality (|2.54p holds for x e {x2,xi), it then also 
follows for x < Xq. 

The inequality (|2.55|) for x > xq follows in the same way. □ 



It follows from Lemma 12.51 that, up to exponentially small corrections, the 
jump matrices on the real line in the RH problem for T take the following form 



/I 

gn(>'2, + -A2.-) Q gri(A4,+ -A2,-) 

10 

VO e"(^^.+ ~^*'-V 



on (—00, X2) (2.56) 
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/gn(Ai,+ -Ai._) Q gn(A3,+ -Ai,_) q\ 

10 

e"(^3+-^3.-) 

\ 1/ 



on (xi, oo), 



(2.57) 



and 



h 



on {x2, Xi). 



The matrices (j2.56p and (|2.57p are in a standard form. By standard ar- 
guments we have that the non-trivial diagonal entries of (|2.56p and (|2.57p are 
rapidly oscillating on the intervals [0.2,(^2] and [ai,/3i], respectively, and they 
are equal to 1 outside these intervals. 

From part (a) of Lemma [TSl we can further deduce that ()2.56p reduces to 



/I I 

e"(^2,4 

I 

Vo ( 

and (|2.57p reduces to 





-) 
1 






1 



•l(A4,+ -A4,_) 



[a2,/32 



(2.58) 



n(Ai 



-) 
1 






o"(-^3,- 





-) 

1/ 



on 



(2.59) 



Outside these intervals, the jump matrices (|2.56p and (|2.57p are exponen- 
tially close to the identity. 



2.4.4 Summary 

Summarizing, we have shown that we can take Fi and F2 so that all jumps 
in the RH problem for T tend to the identity matrix as 71 —> 00, except for 
the jump matrices on the two intervals [aj^(3j\^ j — 1,2, which are given by 
(|2.58p - (|2.59p plus an exponentially small term. 

Observe that each of the asymptotic jump matrices above has the sparsity 
pattern 








X 







X 





X 


X 





X 







X 








(2.60) 



It follows that the 4x4 RH problem for T — T{z) asymptotically decouples 
into two RH problems of size 2x2, one involving rows and columns 1 and 3 
and the other involving rows and columns 2 and 4. The coupling between the 
two RH problems is exponentially small as n ^ 00. Thus we are dealing now 
essentially with two RH problems of size 2x2. The remaining steps in the 
Deift-Zhou steepest descent analysis can then be done in a standard way on 
these decoupled 2x2 problems. This will be described in the next subsection. 
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2.5 Remaining steps of the steepest descent analysis 

2.5.1 Third transformation: Opening of the lenses 

The next step in the steepest descent analysis is to open lenses around the 
intervals {aj,(3j), j = 1,2. This operation serves to transform the oscillating 
jumps on these intervals into exponentially decaying or constant ones. This 
can be done in the usual way [TH] on each of the decoupled 2x2 problems. 
We obtain in this way a new matrix function S = S{z) obtained from T{z) by 
multiplication on the right with a suitable transformation matrix, the precise 
form of which depends on the different regions of the complex plane. 

Let us describe this in detail. First consider the interval [ai,/3i]. We have 
here the matrix factorization 



3Ti(Ai, + -Ai,_ ) 



,ri(Ai,+ -A3, + ) 



1 

,ri(Ai__-A3__) 



1 

g"(A3,+ -A3,-) 

1 

g™(Ai 





,--A3,-) 

1 

-1 



□ "(Ai 



We can then open a lens around [ai, and define 
/ 1 



S{z) 



T{z) 



T{z) 



1 0\ 

10 

_gn(Ai-A3) 1 

V 1/ 
/I 0\ 

10 

2n(Ai-A3) 1 

V 1/ 



in upper lens region around [ai,/3i] 
in lower lens region around [ai, 

(2.61) 



Similarly we open a lens around [a2,/32] and define 

in upper lens region around [a2 , 132] 



S{z) = < 



T{z) 



T{z) 



/I 






Vo 

/I 




vo 



0\ 

1 
1 

gn(A2-A4) 1/ 

0\ 





1 

1 

grx(A2-A4) 1/ 



in lower lens region around [a2, f32\- 



We also set 



S{z) = T{z), outside the lenses. 
The matrix function S = S{z) satisfies the following RH problem 

(1) S{z) is analytic in C \ (R U Ti U r2 U lips of the lenses); 

(2) S satisfies the following jumps: 



(2.62) 
(2.63) 



33 



S+{x) = S-{x) 
S'+(x) = S-{x) 
S+{x) = S-{x) 

S+{x) = S-{x) 



I 










— 




-A2) 


1 




— 1 









U 




u 


/I 




2,+ - 
























Vo 




-1 




/ 


1 













1 







i(Ai-A3) 





1 


V 











/I 













1 















1 


VO 


gn(A2- 


-A4) 






1 0\ 

gn(A3,+ -A2) Q 


1/ 

1 

1 
/ 



on (ai,/3i) 



on (q!2,/?2) 



on the lips of the lens around 



on the hps of the lens around 



["2,^2]- 

On the remaining contours Fi, (—00, a2), (/32, oii) and (/3i, 00), the 
jumps for S are exactly the same as those for T . 



(3) As z ^ 00, we have that 



S{z) = I + 0{\/z). 



As already mentioned before, the entries in positions (1,2), (1,4), (2,1) and 
(2, 3) in the jump matrices for S are all exponentially small when n ^ 00 (use 
parts (a) and (b) of Lemma [T5|) . 



2.5.2 Model RH problem: Parametrix away from the branch points 

We consider now the model RH problem obtained from the RH problem for 
S = S{z) by ignoring all exponentially small entries in the jump matrices. The 
model RH will be defined in the region 

C\([ai,/3i]U[a2,/32]). 

The solution p(°°^(z) to this model RH problem can be constructed in the 
usual way pL6j for each of the two 2x2 problems into which the 4x4 problem 
decouples. This leads to the parametrix 



P(°°)(z) = 



/ 7l(z)+7i ^(z) Q 7l(z)-7i ^(z) Q 

r, 72(z)+72~^(^) n 72(z)-72~^(^) 

" 2 " 2i 

7l(^)-7i ^(g) n 7l(z)+7i ^(z) n 

22 2 

72(z)-72'^(^) n 72(z)+72~^(^) 

2i " 2 







where 



71(2;) := 



z-(3i 
z — ai 



1/4 



, 72(2) := 



z-(32 

Z - Q!2 



1/4 



(2.64) 
(2.65) 



and where we choose the principal branches of the 1/4 powers. 
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2.5.3 Local parametrices around the branch points 



Consider disks around the branch points aj, j = 1,2 with sufhciently small 
radius. Inside these disks one can construct local parametrices p(^"^y) (z) to the 
RH problem for S{z) in terms of Airy functions. Once again, these parametrices 
can be constructed in the usual way [IS] on each of the two 2x2 RH problems. 
We omit further details. 

2.5.4 Fourth transformation and completion of the steepest descent 
analysis 

Define a final matrix-valued function R{z) by 



From the construction of the parametrices it then follows that R satisfies a 
RH problem 

(1) R{z) is analytic in C\T,fi where S/j consists of the real line R, the contours 
Fi and the lips of the lenses outside the disks, and the boundaries of 
the disks around ai. Pi, a^-, and /32, 

(2) R has jumps i?+ = R-Jr on E^j, that satisfy 

Jr{z) = I + 0{l/n), on the boundaries of the disks, 
Jr{z) =1 + 0(e"™l^l), on the other parts of T^r, 

for some constant c > 0. 

(3) R{z) = 1 + 0{l/z) as z oo. 

Then, as in [THl [T71 |TB] we may conclude that 



as rt — > oo, uniformly for z in the complex plane. This completes the RH 
steepest descent analysis. 

2.6 Proof of Theorem 11.31 in the case of large separation 

We will now establish Theorem 11.31 in the case of large separation. The idea is 
that from the asymptotic decoupling of the 4x4 RH problem for large n, there 
follows a similar decoupling for the kernel in (11.23P and hence for the associated 
non-intersecting Brownian particles. 

First assume that x,y £ (aJ,/3*) and consider the correlation kernel (|1.23p 
Kni,n2;ni.n2i which for short we denote by Kn- By virtue of (12. lip we obtain 




5(^)(p(Airy))-l(^)^ 



in the disks around ai, /3i, a2, /?2 
elsewhere. 



(2.66) 




(2.67) 



2'Ki{x — y) 



1 



(0 W2Ay) 0)X^\y)X+ix) 



Wi,2{x) 





V / 
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Using (p:^ and 1^:^ this becomes 



2'Ki{x — y) 
From p.6ip we get 



(O e"^3.+ fe) 0) T+i(y)T+(x) 



/g-nAi, + (x)\ 
g-nA2(2;) 



V y 



K„ix,y) 



(_e"^i.+to) e"^3,+ fe) 0) S'+^(y)5+(x) 



-nAi, + (a;) 
—n\2 (x) 
-nA3, + (x) 



2m{x — y) 

V / 

(2.68) 

Now it follows from (|2.67p by standard arguments (e.g. ff, Section 9]) that 

S^^{y)S{x) = I + 0{x ~ y), as y ^ x 
uniformly in n. Inserting this into (|2.68p yields 

_ n(Ai,+ (iy)-Ai.+ (a:)) , n{\a_+{y)-\3, + ix)) 

ZTn(x — y) 

+ (9(e»(^i,+ fe)-^2(a=))) ^_(9(gn(A3,+ to)-A2(x)))^ y -> 2;, (2.69) 



uniformly in n. Now in the limit when y — > a;, the last two terms in p.69p 
become exponentially small by virtue of Lemma l2.5f b). Then by letting y ~* x 
and using I'Hopital's rule and find 

if„(,,,) = !^iikL:Li^ + o(i) 



■^iPi-x){x-ai) + 0{l) 



27rt(l - t) 

as n ^ oo, where the last step follows from ()2.35p . We conclude that 



lim —Kn{x,x) — - 
n— >oo n 27rt(l — t) 



l—^^{Pl-x){x-al), xeialPl). 



For X e (q!2: P2) obtain in a similar way that ■^Kn{x, x) tends to the semi- 
circle density on the interval [a2,/92]- Thus by (|1.23p we have completed the 
proof of Theorem II. 31 in the case of large separation (recall that T ~ 1). □ 



3 Steepest descent analysis in the case of critical 
separation 

In this section we study the non-intersecting Brownian motions in case of critical 
separation between the endpoints. We will work under the assumption of the 
double scaling regime p.25p - p.26p . This means that we take the endpoints 
ttj ,hj, j = 1,2 fixed such that 

(ai-a2)(&i-62) - (v^+ Vp^)', (3.1) 
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and that we consider the temperature T to be varying with n as 

T = Tn = l + Ln^^'^. (3.2) 

Recall also the assumptions ()1.5p - ()1.8p . 

Before proceeding further, let us recall the main objects needed for the steep- 
est descent analysis. The points j = 1,2 are as defined in (|2.8p - (|2.9p . 
but now with T in (|3.2p not identically equal to 1: 

aj = aj{t) = (1 - t)aj + thj - 4pjTt{l - t), (3.3) 
P] = = (1 - 0%- + + ^ApjTt{l~t), (3.4) 
for j = 1,2. The limiting values for n oo of these points are denoted by (cf. 



a* = a* (0 - (1 - t)aj + tbj - ^ Ap*t{l - t), (3.5) 

p* = p* (t) = (1 - t)aj + tb, + y^4p*i(l - t), (3.6) 

j = 1,2. Recall also the A-functions and ^-functions which are given by ()2.30p 
and (12.35|) . We will denote the limiting values for n ^ oo of these functions by 
X%{z), Cfe(2); fc = 1, . . . ,4. For example, the functions £,'^{z) are given by 



aw 


1 

~ 2t(l-i) 


(- 


(1 


- t)ai - 




+ ((^ 


-al){z 






_ 1 
2t(l-t) 


(- 


(1 


- t)a2 - 


^tb2 


+ ((^ 








1 

~ 2t(l-t) 


(- 


(1 


- t)ai - 


^tbl 


-((^ 


-~al){z 






1 

2t(l-t) 


(- 


-(1 


- t)a2 - 


f t&2 


-((^ 


~al){z 


- P*2)y^' 



(3.7) 



Throughout this section we will assume again that the hypothesis (12. 7p holds, 
but now with strict inequality, i.e., 

0<t< <„it = 7 ^— ^ rr- (3.8) 

(ai - 02) + (bi - 62) 



The time tcrit is now the time where the two ellipses in Figure 2(c) are tangent 
to each other. The case where tent < t < 1 can be handled in a similar way; cf. 
Remark I2. II The behavior at the tangent time icrit itself leads to a multi-critical 
situation which is outside the scope of this paper. 

3.1 Modifications in the choice of xq, Ti and r2 



The main technical tool that made the analysis in Section 12.41 work was the 
existence of a point xq e {^2,0^1) such that the string of inequalities (|2.36p 
holds, cf. Lemma 12.31 It will turn out that in the present context, this can be 
achieved by the point Xq given by the explicit formula 

^l+ft^ (3.9) 



Note that we already encountered the analogue of the point ()3.9p in the proof of 
Lemma 12.31 cf. (j2.37p . In particular, in the proof of Lemma 12.31 we derived in- 
equalities (|2.39p and (|2.40p in case of large separation of the endpoints. Running 
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through the proof of these inequahties, we find that in case of critical separation 
these inequalities become equalities: 



Pi 



and 



{xl-al){xl-l3;) 



P2 



{{l-t){ai-a2)-t{bi-b2)f 



These equalities can be restated in the form 
^[al-xim-xl) = ^J^ 



{{l-t){ai-a2)-t{bi^b2)f 



^{xl~al){xl-l3*) 



'Pi + VP2 



{{I - t){ai - a2) - t{bi - b2)) (3.10) 



{{1 - t){ai - a2) - t(bi - b2)) . (3.11) 



The positivity of the right-hand sides of p.lOp - p.lip follows from our assump- 
tion (IH]) . 

In what follows we will also need the identities 



Oil + fi*2 



IP\ + \JP2 



((l-t)(ai-a2)+i(&i-&2)), (3.12) 
((l-t)(ai-a2)+t(6i-62)). (3.13) 



These identities follow immediately from the definitions p.9p and p.5p - p.6p . 

Lemma 3.1. Under the double scaling regime p.ip - p.2p . the point Xq defined 
in p.9p satisfies the inequalities 



(3.14) 



PROOF. From the definitions p.7p and keeping track of the correct branches 
of the 1/2 powers, we find 



a(2:5)-er(a;S) 



-iil-t){ai-a2)-t{bi~b2) 



2t{l - t) 

-^{x*,^ a2)(xl ~p2) + ^ {xl - ai)(xS - /3i)). (3.15) 
By (|3J0l) and ([3lT|l . this is 

1 ^ ^ 



2t{l - t) 



((l-t)(ai-a2)-t(6i-62))(l 



1 



2t{l - t) 
Similarly, e2*(x5)- a (4) >0 



{{I - t){ai - a2) - t{bi - b2)) 



Pi VPi 



VP2 
> 0. 



□ 
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Note that the relations (|3.14p correspond to the two outermost inequahties 
in (j2.36p . evaluated asymptotically for rt — > oo. By continuity, these inequalities 
then also hold for n finite but sufficiently large. 

In a similar way we would like to have the middle inequality in ()2.36p . i.e., 
^3(a;o) > CI (2^0)- However, this is not the case, since instead we have equality 
^4(2^0) ~ ^3(2^0) = 0- We need a more detailed statement of the zero behavior. 

Lemma 3.2. We have 

\liz)-Xl{z) = c ^'-^f\ oiiz-x*r), z^x*, (3.16) 
where c is the positive constant given by 



((l-t){a,-a2)-t{h-b2))-\ (3.17) 
VP1P2 

PROOF. Recall from Section [5] that we choose the constant k such that 

Xiix*) = X^ix*) (3.18) 

and hence in particular 

•^4(3^0) — ^3(^0)- 

Now we consider the subsequent derivatives of A| — A3 at Xq . We start with 
the first derivative CI (2^0) ~ Ql^^o)- Recall that in the proof of Lemma we 
showed that this derivative is negative provided there is large separation between 
the endpoints. The same proof shows that in case of critical separation, this 
derivative is zero. This means that 

^4(2^0) — ^3(^0) I 

as we already alluded to in the paragraph before the statement of Lemma 13.21 
Consider then the second derivative 

(CD'K) - (QYi^o) (3-19) 

of A4 — A3 at Xq . It follows immediately from p.7p that 

Inserting (|3.10p - (|3.13p . and keeping track of the branches of the 1/2 powers, we 
see that both terms in p.l9p cancel each other when z = Xq and hence p.lPp 
vanishes. 

Consider then the third derivative 



(C|)"K)-(C3)"K) (3.20) 



of A4 - A3 at Xq . We have 



-iiz-al)iz-Pl))-'/') 
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or equivalently 



((^-^)'-(^-"2)(^-/32)) 



which can be further simphfied to 

By inserting (|3.5p - p.6|) this yields 

(a)"(z) = 2p^((z-«D(z-/3r))-3/2 

Using ((m]l) - ([TTT|) . we then find 

l^-i^ + -j=^ + v^)' ((1 - t){a, - a,) - t{h, - h))- 



= 2 



This proves the lemma. □ 

From Lemma 15^ we see that the behavior of Re {Xl{z) — Xl{z)) in the neigh- 
borhood of Xq depends on the different sectors in the complex plane: in partic- 
ular this real part is negative in the sectors 

-| < arg(z-x;5) < ^ < &Tg{z - x*q) < ^ 

and it is positive in the sectors 

-y < arg(z-xo) < --, - < arg(z-xo) < y. 

We will then choose Fi so that it lies in the sectors with negative real part and 
F2 so that it lies in the sectors with positive real part. In particular, note that 
Fi and F2 now both pass through Xq. Just as in Section [21 we will choose these 
curves independent of n. 
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Figure 8: The jump matrices in the RH problem for T = T{z) in the case of 
critical separation. Compare with Figure El 



3.2 Transformations of the RH problem 

We are now ready to describe the steepest descent analysis in case of criti- 
cal separation, assuming the double scaling limit (|3.ip - (|3.2p with n fixed and 
sufhciently large. 

3.2.1 First transformation: Gaussian elimination in the jump matrix 

The first transformation Y ^ X oi the steepest descent analysis is a Gaussian 
elimination in the jump matrix. This is done in the same way as in Section [2. 2) 
except that we choose the contours Fi and r2 in a different way. As explained 
at the end of the previous subsection, these curves should be chosen fixed (in- 
dependent of n) and such that they pass through the point Xq in certain sectors 
of the complex plane. 

3.2.2 Second transformation: Normalization at infinity 

The transformation X i— > T is the same as in Section r2.3l The new matrix valued 
function T — T{z) is normalized at infinity in the sense that T{z) = / + 0(l/z) 
as z — > cx). The jump matrices in the RH problem for T{z) are shown in Figure 

m 

3.2.3 Asymptotic behavior of the jump matrices 

By Lemmas 13.11 and 13.21 and the choice of the contours Fi and F2 , we see that 
the conclusions of Section 12.41 all remain valid for n sufficiently large, except 
that in a neighborhood of the point Xq (through which now both Fi and F2 
pass), the jump matrices are not exponentially close to the identity matrix. 
This means that, except for a small neighborhood of Xq, the RH problem again 
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Figure 9: After ignoring the exponentially small entries of the jump matrices in 
Figure [HI the jump matrices near Xq in the RH problem for T{z) have a 2 x 2 
structure as shown in the figure. 



asymptotically decouples into two 2x2 problems involving rows and columns 
1, 3 and 2, 4, respectively, in exactly the same way as in (|2.58p - (|2.59p . 

The only place where the above decoupling is not valid is in a neighborhood 
of Xq. In fact, since Xq is away from the intervals [ai,/9i] and [0:2, /92] for n 
sufficiently large, the jump matrices on the real line close to Xq are all exponen- 
tially close to the identity matrix. Ignoring these jump matrices, the only jump 
conditions that remain are those on the curves Fi and F2. From Figure [8] we 
see that the latter constitute essentially a 2 x 2 RH problem involving rows and 
columns 3 and 4 only. This leads to the jump matrices shown in Figure [51 

3.2.4 Third transformation: Opening of the lenses 

Around the intervals [ai,/?i] and [a2,/?2]j we are in the region where the RH 
problem decouples into two 2x2 problems. We can then define a transformation 
T I— > S' by opening lenses around these intervals. This is done in exactly the 
same way as in Section 12.5.11 Since these operations have no influence on the 
behavior around Xq (which is our main point of interest here), we omit a detailed 
description. 

3.2.5 Model RH problem: Parametrix away from the special points 

We describe now the solution to the model RH problem where we ignore all 
exponentially small entries in the jump matrices, and where we stay away from 
the special points ai , /?i , 012, /32 and Xq. Since we are considering the region away 
from the special point Xq, we are essentially dealing with two 2x2 matrix valued 
RH problems. The construction of the parametrix P^°°^z) is then exactly the 
same as in Section r2. 5. 21 
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3.2.6 Local parametrix around ai, 0:2, /?2 

In small disks around the endpoints of the intervals [q;i,/?i] and [a2,/32], one 
can construct local parametrices p(^"'y) (2;) to the RH problem in terms of Airy 
functions. The construction is exactly the same as in Section [2.5.31 

3.3 Local parametrix around Xq 

3.3.1 Model RH problem associated to the Painleve II equation 

Our next task is to build a local parametrix in a neighborhood of the special 
point Xq. This will be the main technical part of the analysis. In this subsection 
we will first recall the model RH problem associated with the Hastings-McLeod 
solution of the Painleve II equation. We describe this RH problem and a few 
basic deformations of it, serving to bring it closer to the form we will need. 

We start with the RH problem for the 'vp-functions' associated to a general 
solution of the Painleve II equation [501 HI] ■ To this end, consider a matrix 
function depending on two variables 5* = ^'(C, s), of which the variable s is 
considered fixed at this moment. The jump conditions of 5* = 5'(C,s) in the 
C-plane are shown in Figure 1101 The jumps are determined by three numbers 
si, 82,83 S C which are called Stokes multipliers; these may be any complex 
numbers satisfying the relation 

si + S2 + S3 + S1S2S3 = 0. (3-21) 

The 2x2 matrix valued function 5* then satisfies the following RH problem: 

(1) ^{C, s) is analytic for C e C \ ULoiargC = 7^/6 + kir/S}; 

(2) For e Ufc=o{^^g^ = 7r/6 + kn/i}, ^'(^, s) has jumps as shown in Figure 

(3) As C ^ c)o we have that 

*(C, s)e*(^f'+^f)"^ = / + 0{C^), (3.22) 
where (T3 :— diag(l,— 1) denotes the third Pauli matrix. 

The corresponding Painleve II function is recovered from the RH problem by 
the formula ^ [21] 

q{s) = lim 2iC^'i2(C,s)e'*^C"+<)'^3^ 

— >oo 

The Hastings-McLeod solution to the Painleve II equation corresponds to 
the special choice of Stokes multipliers 

Si = 1, S2 = 0, S3 -1. 

Since S2 = 0, we see that the jump on the imaginary axis in Figure fTUl disappears. 
The jump matrices of the RH problem then reduce to the situation in Figure 
im Note that we have reversed the orientation of the two leftmost rays. 
Now we apply a rotation over 90 degrees, i.e., we set 

$(C,s) s). (3.23) 

Then 4' satisfies the following RH problem: 
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Figure 10: The jump matrices in the ^-plane in the RH problem for = ^{C, s) 
associated with a general solution of the Painleve II equation. The Stokes 
multipliers si, S2, S3 satisfy si + S2 + S3 + S1S2S3 — 0. 



(1) *(C, s) is analytic for C e C \ ULiiargC = ±fc7r/3}; 

(2) On the rays argC = ±7r/3, 4'(C, s) has jump matrix ^ _^ 

On the rays argC = ±2tt/3, vI/(C,s) has jump matrix ^| , see Fig- 
ure [H 

(3) As C ^ 00 we have that 

$(C, s)e(-^«'+"'^)"^' = / + 0(C~'). 

Now we apply one final modification. Define a new matrix function 

M(C, s) := a3$(C, s)e(-^«'+<)"^a3. (3.24) 
Then M satisfies the following RH problem: 

(1) M{C, s) is analytic for ( e C \ ULiiargC = ±fc7r/3}; 

(2) On the rays arg^ — ±7r/3, M{(, s) has jump matrix ( ^ 

On the rays arg^ — ±27r/3, M{(, s) has jump matrix ^ _(^s^3_2gQ , 
see Figure HH 

(3) As C ^ 00 we have that 

M{C,s) = I + 0{C^). (3.25) 
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Figure 11: The jump matrices in the ^-plane in the RH problem for ^ = ^{C, s) 
associated with the Hastings-McLeod solution of the Painleve II equation. The 
Stokes multipliers are si = 1, S2 = 0, S3 = — 1. 



Note in particular that the RH problem for M — M((, s) is normalized at 
infinity in the sense of (|3.25p . Moreover, the O-terni in (|3.25p holds uniformly 
for s in compact subsets of C\7^ with V the set of poles of the Hastings-McLeod 
solution to the Painleve II equation. Of particular interest for us is the fact |2S] 
that P n M = 0, i.e., there are no poles of the Hastings-Mcleod solution on the 
real line. The RH problem for M = M(C, s) will be used in the construction of 
a local parametrix near x^. 

3.3.2 Construction of local parametrix 

Now we construct a local parametrix near the special point Xq . It turns out that 
the construction can be done in a similar way as done by Claeys and Kuijlaars 
in [TT]; see also [12]. We also note that the construction of a local parametrix 
with ^i-functions associated with Painleve II played a role [5] and [6] . 
We are going to construct the local parametrix in the neighborhood 

Us := {z e C| \z-x*a\ < 6} (3.26) 

of Xq, where the radius (5 > is fixed but sufficiently small. 

Recall that for n sufficiently large the jumps near Xq are essentially those 
of a 2 X 2 matrix- valued RH problem involving rows and columns 3, 4, with 
jump matrices shown in Figure [51 To construct the local parametrix near Xq, 
we construct similar jumps via the model RH problem for M(C, s) in Figure [T51 
To this end, we propose a local parametrix P'-^o)(^z) of the form 

^^^'^(^)-^^^^(^)(o M(nV3/(,),,./3,„(,)))> ^^Us (3.27) 
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where p(°°) (z) is the solution to the model RH problem in Section 13.2.51 and 
where we set 

3 



fiz):^ l^-^iXl-X;)iz)j , (3.28) 

We will show in the next lemma that f{z) is a conformal map mapping the 
neighborhood Us of (provided S is small enough) onto a neighborhood of the 
origin in the C-plane. It follows that Sniz) in p.29p is a well-defined analytic 
function in Us-, due to pole-zero cancelation at Xq. Indeed, this follows from the 
fact that both terms in the numerator of p.29p vanish at z = Xq, cf. p.l8p . 

Lemma 3.3. ( Conformal map) The function f{z) is analytic in a neighborhood 
Us of Xq and satisfies 

/(^) = ,w TTT + Oiz - x*r, (3.30) 

2K({1 - t)(ai - a2) - t(bi ~ 62)) 

as z — + Xq, where the constant K is defined in p.49p . 

PROOF. This follows from (|3.28p and Lemma IX^ □ 

From the above definitions (|3?27l) - (|3?29l) we see that if C = n^'^fiz) and 
s — n'^^^Sn{z), then the exponent occurring in the jump matrices in Figure [T51 
reduces to 

^C' - 2sC = n (lf{z) - 2s,,{z)f{z)\ = n{X, - X,){z). 
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Figure 13: The figure shows the jump matrices in the ^-plane in the RH problem 
for M = M{C,s). 



This shows that M(n^/^/(z), n^/^s„(z)) has precisely the required jumps in the 
RH problem near a;Q, cf. Figure [H provided the contours Fi and F2 near a;Q 
are chosen in such a way that they are mapped by the conformal map / to the 
straight lines in Figure [T^ We have indeed the freedom to choose Fi and r2 in 
that way near Xq. 

A technical issue that remains is showing that the RH problem for 
M{n^^^f{z),n'^^^Sniz)) is solvable. This is equivalent with the fact that n^/'^s„(z) 
stays away from V, the set of poles of the Hastings-McLeod solution to the 
Painleve H equation. This will indeed follow from the next lemma. Recall that 
we are working under the double scaling assumption 

T = 1 + Ln"2/3^ 
Lemma 3.4. (Asymptotics of s„{z)) We have as n —f 00, 

n'/'sniz) = , [Fiiz] - F2{z)] + 0{n-'/') (3.31) 

8^/t{l - t)f(z) 

where 

F,{z) = ~ a*) r [{y - a]){y - /?*)] '^'^ dy, j = 1, 2, (3.32) 

Jx* 

and the O-term in p.3ip is uniform for z £ Us- 

PROOF. Since T = 1 + Ln^^^^ and Pj ^ P*j + 0{n^^) it follows from p.3p - 
([^ that 

aj = a* - L^p*t{l - t)n-'^'^ + 0(n"^), (3.33) 
[3j = (3* + L^p*t{\ - t)n-^/^ + O(n-i), (3.34) 
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as n — + cxD. Then by (j3.33p and p.34p we have uniformly for y (z Us, 



-3 ''^^^ 



L^p*t(l-t){P*-a*) 



'2/3 + 0(n-i). (3.35) 



2[{y-a*){y-fi*)]"' 
Since for z £ Us, 

(A4-A3)(z)- r{(4~(3){y)dy, 

Jx* 

ixi-xi){z)= r{a-C3{y))dy, 

Jx* 

we obtain from (|2.35p that 
(A4-A3)(z)-(A:-A*)(z) 

\{y - a2){y - p2)f' - [{y - a;){y - dy (3.36) 



2t{l - t) 
1 



2i(l - 1) J,, 
which by (f?:^ and leads to 

(A4-A3)(z)-(A:-A*)(z) = -— [^^i(^) - F^iz)] n-'/^+0{n-\z-xl)) 

and the O term is uniform for z G Us- 

The lemma now follows because of the definition of s„(z) in (|3.29p and the 
fact that f{z) has a simple zero at z = Xq. □ 

Corollary 3.5. (a) We have and n — > oo , 

) + 0(n-i/3) (3.37) 
where the real constant s is defined in (11.50p . 

(b) //J is sufficiently small, then there is a compact subset IC ofC\V, where 
V denotes the set of poles of the Hastings-McLeod solution of Painleve II, 
such that 

n^/^Sn{z) e IC 

for every z G Us and for all large enough n. 
PROOF, (a) From (|3.3ip it follows that 
n^'s^iz) = J lim^ F,iz) -F,{z) ^ _ ^ ^ /3 
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By (jX^ . (IS3])-(1321), and (PTII)) - (P11|) we have for j = 1, 2, 



and by p.30p wc have 



hm = 2i^((l - t)(ai - 02) - t(6i - 62)). (3.40) 

Multiplying ([STSO]) and ([QO]) we find 

L ^i(^) ; f^(^) = f V^)' ^ .s (3.41) 



8^i(l-t) f{z) 

where for the last equality we used the definition of K, see ()1.49p . and the 
relation pTSOl between L and s. Now (f3^ follows from ([3^ and (|OT|l . 

(b) The equation (|3.37p implies that there exist constants Ki and K2, inde- 
pendent of n and z, such that for z € Us and n large enough, 

\n^/hn{z) - .s\ < Ki\z -~x*o\+ K2n-^/\ 

Since s is real, and since the Hastings-McLeod solution has no poles on the real 
line, part (b) follows. □ 

Now we can check that the local parametrix p(^o) defined by (|3.27p satisfies 
the following 'local RH problem' 



I2 

* 



(1) P(^S) is analytic in Us \ {Ti U Ta); 

(2) On Fi and Fa, P'-^o) has jumps with jump matrices of the form 
where the 2x2 matrices indicated by * are given in Figure [3 

(3) As n ^ 00, 

= p(°°)(z)(/ + 0(n-i/3)), (3.42) 
uniformly for z on the circle dUs. 

The matching condition (3) follows from the definition (|3.27p and the fact that 
M{(, s) = I + 0(1/C) as C ~* 00 uniformly for all s of the form s = n^/^s„(z). 

3.4 Fourth transformation and completion of the proof of 
Theorem D 

3.4.1 Fourth transformation 

Using the global parametrix p(°°) and the local parametrices P^^'^'y) and p(^o)^ 
we define the final transformation S* 1— > P by 

S'(z)(p(^o))~i(2;)^ inside the disk Us around x^, 
R{z)^{ S'(z)(p(A''^y))-i(z), in the disks around ai,/3i,a2,/32, (3.43) 
S'(z)(p(°°))-i(z), elsewhere. 
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From the constructions in this section it foUows that R satisfies the foUowing 
RH problem 

(1) R{z) is analytic in C\Si^ where 'Sn consists of the real line M, the contours 
Fi and F2 outside Us, the lips of the lenses outside the disks, and the 
boundaries of the disks around ai, a2, P2 and Xq; 

(2) R has jumps R+ = R-Jr on E/j, that satisfy 

Jr{z) = I + 0{n^^^^), on the boundary of Us, 

Jr{z) — I + 0{l/n), on the boundaries of the other disks, 

Jr{z) =1 + 0(e"™l^l), on the other parts of S^, 

for some constant c > 0. 

(3) R{z) = 1 + 0{l/z) as z ^ 00. 
Then, as before we may conclude that 



as n 00, uniformly for z in the complex plane. This completes the RH 
steepest descent analysis. 

We will explicitly compute this 0{n^^^^) term in Section[31 

3.4.2 Proof of Theorem 11.31 in the critical case 

We have now completed the steepest descent analysis of the RH problem un- 
der the double scaling regime l|3.ip - l|3.2p . In a completely similar way as in 
Section we can now use the result (I3.44p to show that the non- intersecting 
Brownian particles are asymptotically distributed on the two intervals according 
to two semicircle laws. This completes the proof of Theorem 11.31 



4.1 Formulas for recurrence coefficients 

In this section we investigate the large n asymptotics of the recurrence coef- 
ficients of the multiple Hermite polynomials (Sections 11.71 and II. 8p under the 
assumptions in Section [3l This will lead to the proofs of Theorems 11.91 and 11.121 
We have not yet proved Proposition [TT71 This will be done in Section O but 
anticipating this result, we call the combinations 




(3.44) 



4 Proofs of Theorems 11.91 and 11.121 



1 < i < j < 4 



(4.1) 



the off-diagonal recurrence coefficients, and 



Ci,k+2 



l<hj,k<2, i^j, 



(4.2) 



the diagonal recurrence coefficients. 
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Recall that the Ci^ are the entries of the matrix Yi in the expansion 
Y{z) =(l+^ + o(^\\ diag(z"i , z"^ , z-"i , z-"^) 



as z — !■ oo. 

Following the transformations (|2.13p . ()2.26p . (|2.63p and (|3.43p in the steepest 
descent analysis, we have the following representation for z in the region between 
the contours Fi and F2 



Y{z) = L''T{z)G-'\z)L-'' = L"i?(z)p(°°)(z)G-"(z)L-". 
The following lemma is easy to check. 
Lemma 4.1. We have 

Yi = (Gi + Pf^ + 
where Gi, P^°°^ and Ri are matrices from the expansions as z ^ 00, 



(4.3) 



Gi 



G-"(z)= /+^ + - diag(z"\z"%2""\z 



P^°°\z)=I+^ 



(00) 



i?(z) — + — 



1 



(4.4) 

(4.5) 

(4.6) 
(4.7) 



We are only interested in the combinations (|4.ip and (|4.2p of entries of Yi . 
Since L is a diagonal matrix, the factors L" and L~" in (|4.4p will not play a 
role for these combinations. Also, since Gi is a diagonal matrix (which is clear 
from (j4.5p . since G(z) is diagonal), this does not play a role either. Therefore 
we have for i < j, 



Pi°°^+Ri) (P'T'^+Ri 



and for distinct i,j,k, 



{Pi 



(^) 



P^^^ + i?i 



J, 2 



3-k 



(4.8) 



(4.9) 



In what follows we evaluate P^°°^ and Ri. 
The evaluation of P^*-""-* is straightforward. 

Lemma 4.2. We have 


































V 


-VP2 





/ 



(4.10) 
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PROOF. Recall that p(°°)(z) is defined by means of (PIM)) and ^M^. For 
z — > oo we have the expansions 



Since /3j - = 4:y^pjTt{l - t), see ([2131) -(IS!!), we obtain (|ilUl) from (fTM]) 
and gSl). □ 



4.2 Second term in the expansion of the jump matrix for 

Riz) 

We use the clockwise orientation on the circle dUs around Xq . Thus the outside 
of the circle is the +-side and the inside of the circle is the — side. The matrix 
valued function R{z) in (|3.43p then satisfies the jump condition 



i?+(z)-i?_(z)A(z), 
on dUs, with jump matrix A(z) given by 
A(z) =p(=^S)(2)(p(°°))-i(z) 



z e dUs 



(4.11) 







M{n^/^f{z),n^/^s. 



„(.))) (p'-'r'(^)^ (4-12) 



see p.27p . To prepare for the evaluation of Pi, wc first compute the second 
term in the expansion of the jump matrix A{z) as n — > oo. 

Lemma 4.3. (Asymptotics of jump matrix for R{z) ) The jump matrix A(z) in 
()4.1ip . (|4.12p has the asymptotics 



A(2) = / + A(i)(z)n-i/3 + (9(^-2/3)^ 
as n —> oo, where 



z e dUs, 



(4.13) 



Aa)(.)^^p(-)(.) 



/O 



u{n^l^Sn{z)) 

VO g(n2/35„(^)) 






-g(n2/3s„(z)) 
-u{nV^Sn{z))) 



(p(-))-i(z), 



(4.14) 

uniformly for z on the circle dUg. Here q{s) denotes the Hastings- McLeod so- 
lution to Painleve H, u{s) denotes the Hamiltonian 



(4.15) 



u{s) := {q'{s)f~-sq\s)~q\s), 
and f{z) is the conformal map in p.28p . 

PROOF. First, we calculate the large n asymptotics of the matrix M{n^^^ f{z), 
n^^^Sniz)) which occurs in the expression (|4.12p . Recall from Section [5. 3. II that 
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M(C, s) is defined as an easy modification of the model RH matrix ^(C, s) asso- 
ciated to the Hastings-McLeod solution to the Painleve II equation. It is known 
(see e.g. [5T]) that the asymptotic expansion (|3.22p of '^'{C, s) in powers of 
can be refined to 



2^C -u{s)) VC 



oo, 



where q{s) denotes the Hastings-McLeod solution of Painleve II and u{s) denotes 
the Hamiltonian (|4.15p . Following the sequence of transformations ^' i— > 4' i-^ Af 
in Section [3.3.1l we obtain 

It follows that 

/lff«l/3 f(^\ „2/3„ . ^^ _ r , 1 ( u{n^/^ S,,{z)) ~q{n^/^ Sn{z))\ 1/3 

+ 0(n-2/3)^ (4_^g) 

as n — » oo, uniformly for all z on the circle dUs- 

Then (liTl| follows from (liT^ . and fil^ . □ 

Note that A'^^'(z) depends on n, as is obvious from the appearance of 
n^^^Sniz) in (I4.14p . Also p(°°) depends on n. As n ^ oo, the n-dependent 
entries have limits, and therefore we can still obtain from the expansion (|4.13p 
of A(z) a similar expansion of the RH matrix R{z) 

R{z) = / + {z)n-^/'^ + 0(71^2/3) as n ^ oo. (4.17) 

The coefficient i?(^^(z), in this expansion can be found by inserting (|4.13p 
and (|4.17p into the jump condition R+{z) = R-{z)I^{z) and collecting terms of 
order n^^^'^. This leads to the following additive RH problem for i?'^^: 

(1) i?(i)(z) is analytic in C \ dUs, 

(2) R^^\z) = R^l\z) + A(i)(z) for z G dUs, 

(3) i?(i)(z) = 0{l/z) as z -> oo. 

Note that i?'^^ also depends on n. 

The jump matrix (|4.14p has a simple pole at z = Xq. As in [HI US], the RH 
problem for R^^^ then has the explicit solution 

f ^-;rRes(AW(z),z = xS)- A(i)(z), z£Us 
i?(i)(z)=<^ ^-^0 (4.18) 
jRes(A(i)(z),z = a;S), zeC\Us. 

\ Z 

As in other works (see e.g. PTl f5H|) the expansions of R{z) as z ^ oo and 
n ^ oo commute with each other. It thus follows from ()4.17p and (I4.18P that 

i?i = lim zR^^\z)n-^'^ + 0{n-^l^) 

= Res(A(i)(z), z = x*)n-^/^ + ©(n^^/S). (4.19) 
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4.3 Evaluation of Ri 

In view of (|4.19p the evaluation of i?i comes down to the determination of the 
residue of A^^^ at z — Xq. The resuh is the following. 



Lemma 4.4. We have 
Ri = K 



as n —> oo, 



(4.20) 



where K is the constant defined in (|1.49p and where the 2x2 blocks A, B, C , 
D are given by 



~u{s) -~q{s) 
q{s) u{s) 



A = t{bi - 62) 



C 

D ^ (1 - t){ai - 02) 

with the value s given by (|1.50|) . 
PROOF. From (I3.30p we obtain 

z — x" 



u{s) -q{s) 
q{s) -u{s) 



(4.21) 
(4.22) 
(4.23) 
(4.24) 



lim 



f = X((l-i)(ai-a2)-t(6i-62)), 



z^x* 2f{z) 
and from p.37p we have 

q{n'/hnix*o)) = q{s) + 0{n-'/'), u{n'/'s„{x*)) = u(,s) + 0{n''^') 
as n — > 00, so that (|4.14p yields 

Res(A(i)(z), z = X*) = K{{1 - t){ai - 02) ~ t{bi - 62)) 
/O \ 



xp(-)(xS) 





u{s) -q{s) 
\0 q{s) ~u{s)J 



(p(°°))-i(xS) + 0(n-i/3)_ (4 25) 



To evaluate P'-°°''{xq) we observe that by (|2.64p we need the entries 7j(a;Q)± 
7j(a;g), for j ~ 1,2. We compute the squares of these 

{^,{x*) ± ij'ixDf = 7|(a;S) ± 2 + 77^(4) 
= ±2 



' ^0 Pi 



x^ — a 



^0 ~ /^j 



±2 + (-l)^ 



2x1 - "3 - 



J = 1,2, (4.26) 



where we used (|2.65p . We emphasize that we use ^f- to denote the positive 
square root of a positive real number. Now recall that aj and j3j are varying 
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with n, and satisfy aj — a * + 0(n-2/3)^ = ^* + 0(„-2/3). For the starred 
quantities we can use the identities (|3.10p - p.l3p . so that (|4.26p yields 

(1 - i)(ai - 02) - t(6i - 62) 

Since < 72(2^0) < 1 < 71(2^0) (e^s can easily be seen from (j2.65p and the fact 
that a2 < /32 < < «! < we obtain after simple calculations 

7j(a^5) + 77'(4) _ v/(r-0(Qi-a2) _^ n(„-2/3. u 27) 

2 V(l-^)(ai-a2)-t(6i-62) 

2 V(l-i)(«i-«2)-i(&i-62) 

(4.28) 

for J = 1,2. 

Putting ([42711 and (14281) into (|2J4ll we find 



V(l-t)(ai-a2)-i(foi-&2) 



/V(l-t)(ai-a2) - 62) \ 

V(l -0(ai -02) i^UbT^ 

i^tih - &2) v/(l-0(ai-a2) 



\ -i - 62) v/(l-^)(«i -«2)/ 

+ 0(n-2/3)_ (4,29) 

Using the fact that P^°°\xq) has a 2 x 2 block structure with blocks having 
determinant one, we then also obtain 

(p(-))-i(xS) ^ 



^{I~t)iai-a2)~t{bi-b2) 



/y/{l-t){ai-a2) iy/tibi - 62) \ 

-02) -Wtibi - h) 

~i^t{bi - 62) V(l -0(ai -^2) 

V t^t{bi - 62) v/(l-0(ai -02)/ 

+ 0(n-2/3)_ (4 30) 

Then we insert ()4.29p and (|4.30p into (|4.25p and after simple calculations we 
obtain (|420l)-(|424l), see also (|449l) . □ 



4.4 Proofs of Theorems [TH] and [TT^ 



Having Lemmas 14.21 and 14 . 41 it is now straightforward to compute the asymptotic 
behavior of the recurrence coefficients from (|4.8p and (|4.9p . 

As an example, let us compute the expression in (ll.40p . By ()4.9p and 

from the fact that (^'i°°'')i,j — whenever i + j is odd, see ()4.10p . we have 



'■,3 

Cl,4 {Ri)iA 



(4.31) 
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By (11201) and ^^-^^ we have 

(iil)l.4 Si,2 V^(l-^)(ai-«2)(fcl-fe2) 

V ai-a2z^t(l-t) ^ ^ ^ ^ 

Furthermore, by (li7TU| and 

(-Pi'°°^)2,4 + (i?i)2.4 = ^Vp^VrtO^ + 0{n-^/^), 

which since P2 = P2 + and T = 1 + 0(n^^/'^), reduces to 

(/'l^°"^)2,4 + (i?l)2,4 = iv^yt(r^+ 0(n-l/3). 

Using this and ((13^ in fOT]) . we obtain 



Cl,2C2,4 _ ^ /^, 61-62 , ^^„_i/3 



i\/P2^ i + 0(„-V3), (4.33) 

ci,4 V ai - 02 

which proves the second formula in Theorem II. 121 The other formulas in The- 
orems 11.91 and 11.121 are established in a similar way. □ 

Remark 4.5. (The case where ^crit < t < 1) Recall that the above derivations 
have all been made under the assumption that < t < icrit- The case where 
^crit < t < 1 can be handled by means of similar calculations; see also Re- 
mark 12.11 Alternatively, one can immediately reduce this to the case where 
< t < icrit by virtue of the involution symmetry in Corollary 15.131 After a 
straightforward calculation, one sees then that the expansion in Lemma l4.41 and 
hence the conclusions of Theorems 11.91 and I1.12| remain valid for tcrit < t < 1 
as well. 

5 Proofs of Propositions 11.71 and 11.81 

In this final section, we prove Propositions 11.71 and 11.81 We start by consider- 
ing recurrence relations for general multiple orthogonal polynomials. Next we 
specialize these results to the multiple Hcrmitc case (i.e., the case of Gaussian 
weight functions (I1.18p - (|1.19|) ). An important tool will be the Lax pair satisfied 
by the multiple Hcrmitc polynomials. The compatibility condition of this Lax 
pair allows us to derive more detailed information for the multiple Hermite case. 
For example, we show how this condition implies certain scalar product relations 
between row and column vectors in Yi; these turn out to be essentially given 
by the numbers of Brownian particles associated with the different starting or 
ending points. 

This section is organized as follows. Section 15.11 discusses recurrence rela- 
tions for multiple orthogonal polynomials in a general setting. The rest of the 
section is devoted to the multiple Hermite case. In Section we discuss a dif- 
ferential equation satisfied by multiple Hermite polynomials. Section 15.31 deals 
with the compatibility condition of the Lax pair for multiple Hermite polynomi- 
als. Sections 15 .41 and [5.51 deal with the induced scalar product relations. Finally, 
Section 15.61 completes the proof of Propositions 11.71 and 11.81 
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5.1 Recurrence relations for general multiple orthogonal 
polynomials 

We start by discussing the recurrence relations for general multiple orthogonal 
polynomials, leading in particular to the proof of the expression for the off- 
diagonal recurrence coefficients as given in Proposition II. 71 

It is a classical result that orthogonal polynomials on the real line satisfy a 
three-term recurrence relation. In this subsection, we discuss the p + q+1 term 
recurrence relations for the multiple orthogonal polynomials of Definition 11.41 

As in the classical case p = q — 1 [11], it turns out that the recurrence 
relations can be conveniently derived from the Riemann-Hilbert problem in 
Section 11.41 To this end, we introduce the next two terms in the z ^ oo 
asymptotics of the matrix Y{z) = Y„,m{z) in p.22p : 

Y{z) = (i+Xl + Yi + o(1-]) diag(z"^ , . . . , z-"^ , . . . , z-"'). (5.1) 



Note that the matrices Yi (^1)11.111 and I2 = (5^2)n,m are constant (indepen- 
dent of z). 

Our goal is to find recurrence relations between the RH matrices with multi- 
indices n -I- Gfe. m + e; and n, m, respectively. 

Convention 5.1. Above we use and e/ as row vectors. In some of the matrix 
calculations that follow, we also use these standard basis vectors, but then we 
see them as column vectors, as is usual in matrix algebra. For example e^-e^ 
where is the transpose, denotes a matrix whose only nonzero entry is a one at 
the (j, /c)th position, while eje^ = 6j^k is the scalar product of the two vectors. 
So we introduce the convention that is a row vector when used as a multi- 
index. In matrix computations it is considered as a column vector. We trust 
that this will not lead to any confusion. 

To derive the recurrence relations, assume some fixed k G {l,...,p} and 
I e {1, . . . , g}. From the fact that the jump matrix of Yn,m{z) on the real line 
is independent of n, m (cf. (|1.20p . (|1.2ip ). it follows that the matrix valued 
function 

Un',m{z) Yn+ek,m+e,{z)Y^^^{z) (5.2) 

is entire. From (|5.2p and the normalizations at infinity in ()5.ip . it follows that 
for z —^ oo we have 

Unfmiz) = ('^+ -(5^l)n+efc,m+e, 




Here the middle factor in (|5.3p is just the identity matrix with its fcth and 
(p + l)th diagonal entries replaced by z and 0, respectively. By Liouville's 
theorem it then follows that C/^ m is a polynomial and so by 



U^fmiz) ^ zetel + ejej + (yi)„+e,,m+e,efcej-efeej(ri)„,m. (5.4) 

Summarizing, we obtain the following proposition. 
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Proposition 5.2. (Forward matrix recurrence relations) Let k G {1, . . . ,p\ and 

I G {1, . . . , g}. Then we have the matrix recurrence relation 

^^n+efc,m+e, (z) = t^n.m (^)^n,m (2) , (5.5) 

with U^f^{z) given by (|5.4p . 

Evaluating the different entries of the matrix recurrence relation (|5.5p leads 
to a set of scalar recurrence relations for the MOP. In fact, since the multiplica- 
tion with the matrix C/^'^ in (|5.5p is performed on the left, it easily follows that 
the fcth component functions Ak{x) of the MOP, k = 1, . . . ,p satisfy all the same 
recurrence relations. So the recurrence relations for the individual components 
Ai(x), . . . , Ap{x) can be conveniently stacked into vector recurrence relations 
for the vectors A{x). 

Let us illustrate this for p = q = 2, k^l = l. Then the recurrence relation 
((Oil becomes 



+ei,m+ei 



V 



Cl.l - Ci,i 


-Cl,2 


-Cl,3 


-Cl,4\ 


C2,l 


1 








C3,l 











C4,l 








1 / 



i"n,m, (5.6) 



where we use Cij and Ci,j to denote the entries of (yi)n,m and (ll)n+ei,m+ei, 
respectively. Evaluating the first row of this expression and using (|1.37p leads 
to the vector recurrence relation 

-^n+ei+ei,m+ei ~ V-' "-l,! "-l^^-^n+ei ,m ''l,2-^n+e2,m 

- c^^^{-2i:i)K^l'^_,^ - c,A-27Ti)AZ:^^,^. (5.7) 

This is the desired five-term recurrence relation. Note that the factors — 27rz are 
due to the diagonal matrix D in Theorem II. 61 

Instead of the first row, one could also evaluate one of the other rows of 
(|5.6p . Doing this leads to an additional set of (mostly trivial) relations. 

The relation (|5.7p has the undesirable feature that it involves vectors of 
MOP with different types of normalizations. To express everything in terms of 
one type of normalization, for example the (//, 1) normalization, we need the 
transition numbers between different types of MOP normalizations. 

Definition 5.3. (Transition numbers) Let n, m with |n| = |m| -I- 1 be such that 
the multiple orthogonal polynomials exist. We define the transition number 
tn,rn~*^'^^ G C to be the nonzero constant so that the vector relation 

A _ JII,k^I,l) . „x 

-^n,m — '^n,m ^n,m W-°7 

holds. In a similar way we define the transition numbers tn,in^^^'''\ ir^m^^ ^'', 
and t'i^irT^^'^^ ■ 

The transition numbers are contained in the matrix Yi as follows. 

Proposition 5.4. (Interpretation ofYi via transition numbers) Let n, m with 
|n| = |m| be such that the solvability condition holds. Let D be the diagonal ma- 
trix in Theorem ] 1.6[ Then the off-diagonal entries of the matrix D~^{Yi)n,mD 
can be expressed as transition numbers between different types of normalizations 
of MOP. More precisely, the entries of D~^{Yi)n.mD are given as follows: 
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• If ^ 1^ k,k < p with k ^ k then the (fc, k) entry is t^n+e^m'^\ 

• If 1 < k < p and I < I < q then the {k,p + I) entry is ^n^+et^ '^' 

• If 1 < I < q and 1 < k < p then the (p + I, k) entry is t''^'^^^^'^\ 

• If ^ l£ It^- l£ q with I 7^ I then the {p + l,p + 1) entry is ^n^i^ej'^ ■ 

PROOF. This follows easily from (|5.ip . p.37p . and from the definition of the 
normalizations. The result is most easily seen by evaluating (jS.ip column by 
column; we do not provide the details here. □ 



Let us illustrate Proposition 15.41 for the case where p = q = 2. In this case 
the proposition asserts that 



(yi)„,m = D 



( ^ ,(//,1^7/,2) AII,\-,I,\) {U,\^1.2) \ 

n+ei ,m n+ei ,m n+ei ,m \ 

(11,1^11.1) (11,2^1.1) (11,2^1,2) 

n+e2 ,m n+e2 ,m n+e2 ,m 



D-\ 



(1,1^11,1) (1,1^11,2) 

n,m— ei n,m— ei n.m— ei 

\ ,(/,2-.//,l) (/,2^//,2) ,{I.2^IS) I 

\ n,m— e? n.m— 62 n.m— 62 / 

(5.9) 

where the diagonal entries denoted with * are unspecified, and where the diag- 
onal matrix D = diag(l, 1, — 27ri, — 27ri). 

We use the information in Proposition l5.4l to obtain more meaningful expres- 
sions for the recurrence coefficients. To this end, consider again the five-term 
recurrence relation (|5.7p . This relation yields a connection between several MOP 
with different multi-indices and different types of normalizations. To express ev- 
erything in terms of type (//, 1) normalizations it suffices to multiply with the 
appropriate transition numbers from (|5.9p . Then (|5.7|) transforms into the new 
relation 

-^n+ei+ei,m+ei ^ + '^14 '^14 j-'^n+ei ,m Ci^2C2,l A^^.^^ ,^ 

- Ci^3C3,l A;,,m-ei - ClACi,lA'n,m-e2 > (5-10) 

where we recall that Cij, Cij denote the {i,j)ih entry of (Yi)n,m, (i^i)n+ei,m+ei, 
respectively. Note in particular that the factors — 27rz in ()5.7p have disappeared 
in ([STO]) 

The relation (|5.f Op was derived for the special case where p = q = 2 and 
k = I = 1 in the recurrence relations. In a similar way, one can find the 
recurrence relation for general p, q, k, I. 

Proposition 5.5. (Forward recurrence relations) For any k G {1, . . . ,p} and 

I Cz {I, . . . ,q}, we have the p + q + 1 term recurrence relation 

-Vc ~c^ A^^^'^'> -Vc ^c ~ A^^^'^'> f5 11) 

k^k I 

where we use Ci_j and Ci,j to denote the entries of (11)11. m and (i^i)n+efc.m+e! ; 
respectively. In the above sums, it is assumed that k runs from 1 to p while I 
runs from I to q. 



59 



Note that ProDOsition l5.5l implics the recurrence relations in Proposition ! 1.71 
except for the exphcit form of the first term in the right-hand side of each of 
(jOgil - pri^ . which wiU be estabhshed in Section [OX We note that Proposi- 
tion l5.5l is vahd for general weight functions wi^k and W2,i, not only for Gaussian 
weights. 

Definition 5.6. (Off-diagonal and diagonal recurrence coefficients) We will 
refer to the coefficients of the form Ck,ici^k in (|5.1ip as the off-diagonal recurrence 
coefficients. The coefficient Ck,k — Ck,k hi (|5.1ip will be called the diagonal 
recurrence coefficient. 

We use the terminology 'off-diagonal' and 'diagonal' recurrence coefficients 
in analogy with the case of classical orthogonal polynomials p — q — 1. Indeed, 
in the latter case these recurrence coefficients correspond precisely to the off- 
diagonal and diagonal entries of the Jacobi matrix. 

We collect the off-diagonal recurrence coefficients of Definition 15.61 in an 
upper triangular matrix 



H = 



/o 



Vo 



Cl,2C2,l 





Cl,p+qCp+q^l \ 



'^p+q-l.p+qCp+q.p+q- 

/ 



(5.12) 



The matrix H = (i?)n,m is the Hadamard product (entry-wise product) of the 
strictly upper triangular part of Yi = (Yi)n,m with the transpose of the strictly 
lower triangular part of Yi. For example, when p = q ^ 2 we have 



H = 








Cl,2C2,l 







Cl,3C3,l 
C2,3C3,2 






Cl,4C4,l^ 

C2,4C4,2 

C3,4C4,3 

/ 



(5.13) 



The diagonal recurrence coefficient Ck,k — Ck.k in (|5.1ip is equal to a difference 
of two entries, one from (yi)n+efc.m+ei and the other from (Yi)n,m- It is possible 
however to express this quantity completely in terms of the multi-indices n, m. 
The resulting expression is then somewhat more complicated, involving entries 
of both the matrices Yi — (yi)n,m and Y2 — (l2)n,m- 

Lemma 5.7. (Diagonal recurrence coefficients) The diagonal recurrence coef- 
ficient Ck,k — Ck,k in (|5.1ip can be expressed as 



Ck,k — Ck,k 



(yi)fe,fc + 
(^2)fe,p+; 



{Y2)k,p+i - {yi )k,p+i 



{Yi)k,p+i 



Ek 



k^k ^kl^k,p+l 



E, 



/ ^k,p+l^p+l,p+l 



Ck,p+l 

where in (j5.14p the index k runs from 1 to p while I runs from 1 to q. 



(5.14) 



The proof of (|5.14|) is similar to the proof in the case of classical orthogonal 
polynomials, see [TBI Section 3.2]; we omit the details. A more concise expression 
in the case of Gaussian weight functions will be derived in Section 15.6.11 
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In addition to the forward recurrence relations in Propositions 15.21 and 15.51 
one can also run these recurrence relations in backward order. To this end, 
consider the matrix function 



Unlmi^) '■— ■^n,m(2)^n+efc,m+e, (^)- 

(Note that this is the inverse of the matrix m(^) 
approach above one finds: 



By copying the 



Proposition 5.8. (Backward matrix recurrence relations) Let k € 
and I e {1, . . . , 5}. Then we have 

where 



(5.15) 



f^n;m(^) ^zBp+ie^^i+ ^ e^e^- +(yi)„,mep+;ep+;-ep+,ep+;(yi)„+e,,m+e,. 

(5.16) 

Moreover, the matrix Un'mi^) inverse of the matrix U^''^{z) in ()5.4p . 



Proposition 5.9. (Backward recurrence relations) For any k € {1, . . . ,p} and 

I £ {1, . . . , we have the p -\- q + 1 term recurrence relation 

k 

where Cij,Cij denote the entries of {Yi)n.m and (li)n+efc,m+e, , respectively. 
For example, when p = q = 2, k = l = l the above relations become 








Cl,3 


' \ 





1 


C2,3 





-C3.1 




Z - C3,3 + C3,3 


-C3.4 


V 





C4,3 


1 / 



n+ei,m+ei (-z). 



and 



''^-n.m-ei - + ^3,3 " C3,3jA„+g^^ - C3_iCi_3A„^g^_,_e^^_,_e^ 

- ~ ~ ~ a(^'1) 

- C3,2C2,3An+ei+e2.m+ei " ^3,404,3 A„+g^ . 

5.2 Differential equation for multiple Hermite polynomi- 
als 

In this subsection we derive a differential equation for multiple Hermite poly- 
nomials, i.e., assuming that the weights are Gaussian as in (|1.18p - (|1.19p . The 
differential equation was already described in jl3] . 
For general p and q we consider the weights 



W2,l 



,q, 
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with iV > a certain constant. 

We introduce a modification of the RH matrix Y{z). 

Definition 5.10. Define the matrix fmiction 

vl/(z) := Y{z) diag(/i(z), . . . , (5.18) 
where we used the following functions fj, j — 1, . . . ,p + q: 



/fe(z) = 6^^^^"'"^^^"*^'"='=^ k^l,. 



fp+liz) ^ e-T^'^^ l = l,...,q. 



(5.19) 



The reason for defining the matrix function ^! in Definition 15.101 is that it 
has a constant jump (independent of x) along the real line. More precisely, it 
satisfies the following Riemann-Hilbert problem: 

(1) ^> is analytic on C \ M; 

(2) For a; e M, it holds that 



*+(x) = *_(x) 



/, 



where Ipxg denotes the p x q matrix having all entries equal to one; 
(3) As z — > oo, we have that 

*(z) = (/ + 0(l/z))x 

diag(z"Vi(^), • . • , ^"''/p(^), z-^Vp+il^), . • • , z-"'/p+,(z)). 

Since the function ^' is defined from Y by the multiplication on the right 
with an n, m-independent matrix (recall that N is assumed to be constant), it 
satisfies exactly the same recurrence relations as the original RH matrix Y , cf. 
Section 15.11 Moreover, from the fact that ^' has a constant jump matrix we 
obtain that the matrix valued function 

V„,m(2) := *'„^„(z)*„,™(z)-\ (5.20) 

where the prime denotes the derivative, is analytic in the full complex plane. 
By using ([FTTjl . (|5?T8l) . ((57T9| . and ((5?20| . we obtain that as z ^ oo, 



z 



N 

diag(-z + (1 - t)ai, . . . , -z + (1 - t)ap, -tbi, ... , -tbq)x 



t{l-t) 



1+ (^i)"-'" + ...] +0(l/z). (5.21) 



By Liouville's theorem, it then follows that Vn,m is a polynomial and hence we 
obtain from l|5.2ip that 

N fzlp^Da -Ci2 
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where we define the diagonal matrices 

Da := (1 - t) diag(ai, ... ,ap) , Db := t dia.g{bi, . . . ,bq) , (5.23) 
and where we define C12 and C21 by the partitioning of (Yi)n,m into blocks: 

(ri)„,™^:f^" ^'A, (5.24) 

with the diagonal blocks being of size p x p and q x q, respectively. In other 
words, we have put Cn := [Qj]ij=i^...,p, C12 := [cij]i=i^...,pj=p+i^...,p+q, C21 := 
[cij-]i=p+i^..._p+g and C22 [cij]ij=p+i^...^p+q- 

Multiplying both sides of (j5.20|l on the right with ^n,m(z), we conclude: 

Proposition 5.11. (Differential equation for multiple Hermite polynomials) 
The multiple Hermite polynomials satisfy the matrix differential equation 

n,m 

(z), (5.25) 

with Vn^-a,[z) given by (102]) - (f5?24l) . 



5.3 Lax pair and compatibility conditions 

In this subsection we investigate the compatibility condition between the differ- 
ential equation (|5.25|) in Section [5T2l and the recurrence relations (|5.5|) . ()5.15p in 
Section 15.11 These relations together constitute the Lax pair for multiple Her- 
mite polynomials and the compatibility yields nonlinear difference equations the 
recurrence coefficients. 

Throughout this subsection we will use the partitioning of (li)n,m in (|5.24p 
and we will also use the similar partitioning 

(FO„+e.,^+e, = f^" (5.26) 
\021 022/ 

where again Cn has size p x p, and so on. 

The derivation of the compatibility conditions below will be merely a straight- 
forward calculation. The reader who wishes to avoid this kind of calculations 
could take these results for granted and move directly to Section [53] where we 
discuss the (elegant) scalar product relations induced by these compatibility 
conditions. 



5.3.1 Forward relations 

Fix k E {1, . . . ,p} and I E {1, . . . ,q}. The compatibility condition between the 
differential equation (|5.25p and the forward recurrence relation (|5.5p is 

(t^n,m)'(^) = K.+e.,™+e, (zW^U^) - U^U^W^.M, (5.27) 

where the matrix Vn miz) is given by ()5.22p . Vn+ek m+e, (z) is given analogously 

by 

N f zip -Da -Cxi 
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and U^''^{z) is given by (|5.4p . Using the partitionings (|5.24p and (|5.26p we 
rewrite the latter in block form as 

Un'Li^) = P ^ " l)efeel' + Ciiefce^ - efceJCn -efce^Ci2\ ^ 29) 
V C'2iefee^ / - e;ef / ' 

where the top left block is of size p x p, and so on. In the right-hand side 
of (|5.29p we use our convention that and e; denote column vectors, recall 
Convention 15.11 

Inserting (lO^ . ([S?^ . and ([S?^ into ffd^Tl}, we find 

e^e^ OA N 

OJ t(l-i) 
zl - Da -Ci2\ fl +{z~ l)ekfil + CiiBkel - efceJCn -efce|'Ci2 

(521 -Db / V (5*216^6^ / - e;ef 

/ + (z - l)efee^ + CiiefeeJ - efce^Cii -ek^iCx-A ( zl - Da -Ci 



C2iefce^ I-eief J \ C21 D, 

The left-hand side is independent of z, and therefore all terms on the right- 
hand side involving z or cancel out (this can also be checked by direct calcu- 
lation). Hence the equation reduces to 

i(l - t) fekel 



N \ 

-Da -Ci2\ (I - ekel + Ciiefcel" - efce^Cn -efeeJCi2 
C21 Dt J\ C2iekel I - eief 

'I - ekel + CiiBkel - efceJCn -efee^Ci2\ f-Da -Ci2\ 
C2iekel I-eief / V '^21 Dt J ' 

From the identity (|5.30p we obtain the following for the respective blocks 

• (1, 1) block of (jOni) : 

t{l-t) 



(5.30) 



N 



-efeGj, = [DaJ-eke^.+Cueke,.~ekef,Cii]+Ci2C2ieke^.~ekef.Ci2C2i, 



where [•, •] denotes the usual commutator of square matrices. Since diag- 
onal matrices commute with each other, we get 

^^i^^-^e^e^ = [Da,Ciiekel - efee^Cn] + C\2C2iekel - efce^Ci2C2i. 
This implies for the (fc, k) diagonal entry 

{Ci2C2i)^^^ - (C12C21),,, - (5.31) 

for the {k, k) entry, where we also use ()5.23p . 

fCi2C2i)^ =(l-t)(afc-a^.)(Cii)^ , (5.32) 

V / k.k \ / k.k 
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and similarly for the (fc, k) entry, 

= -(1 - t){ak - a^) {Cn),rk . (5-33) 

for any k £ {1, . . . ,p} with k ^ k. 
(1,2) block of (IOO)1 : 

= -(1 - t)akekelCi2 + Ci2{I - e^eH 

- (/ - ekel + Ciiekel - 6^6^ Cii)Ci2 - ekelCuDt- 

Evaluating the lih column of this equation, one obtains 

= --(l-t)afeefee^Ci2e;+0 

- (/ - GfeGfe + Ciiekel - efcGfe Cii)Ci2e; - tbiekelCi2ei, 

which can be rewritten as 

(/ + {tbi + (1 - t)ak - l)ekel + CnefceJ - efce^Cn) C^ei = 0. 

(5.34) 

Similarly one can evaluate the Ith column to obtain 

(^I + {tbj + (1 - t)ak - l)efee^ + CnGfee^ - efce^Cn j Ci2ej = (7126^ 

(5.35) 

for any / S {1, . . . ,q} with I ^ I. 
(2,1) block of dODH : 

= -C'2i(/ - e^e^ + CiiGfeefc - efce^Cn) 

- DbC2iekel - (1 - t)akC2iekel + (/ - e;ef )C2i. 

Similarly to the case of the (1,2) block entry, one can now evaluate the 
different rows of this equation. This yields 

eTC2i (/ + ((1 - t)ak + tbi - l)ekel + Cnefcef - efcef Cn) = 0, 

(5.36) 

and 



C21 (/ + ((1 + t)ak + tbj - l)ekel + Cueuel - efcel^Cn) = ef C21, 

(5-37) 

for any / G {1, . . . , (7} with I ^l. 
• (2, 2) block of ((Oni) : This is just the trivial relation = 0. 
Note that by comparing (|5.35p and (|5.37p . one finds the relations 

(C2lCi2).^„= (C2lCi2)^ (5.38) 

for any I £ {1, . . . , g} with I ^ I. 
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5.3.2 Backward relations 



The relations in Section [5.3. II constitute in fact only half of the available com- 
patibility relations. The second half is obtained from the differential equation 
(|5.25p and the backward recurrence relation (|5.15p . Alternatively, these relations 
may be obtained from the forward relations by using the involution symmetry 
to be described in Corollarv l5.13l 

We will not derive all these relations in detail. Instead, we list only the 
equations that we will need in the sequel. These are the following analogue of 

(c2iCi2)^ ^ - - fc^, (5.39) 

and the following analogues of (|5.33p and (|5.38p : 

iC2lC,2)li = t{bi - hj) {C22)U , (5.40) 
for any I E {1, . . . ,q} with I 7^ /, and 

(^12(721).. = (Ci2C2i)^,S (5.41) 

for any k e {1, . . . ,p} with k ^ k. 
5.3.3 Involution symmetry 

As an aside, this might be a good place to write down explicitly the involution 
symmetry which we referred to in Section 15.3.21 as well as in Remarks 12.11 and 

EH 

We consider the involution 

We will temporarily denote the RH matrices corresponding to the original and 
reversed RH problem with ywi,n;w2,m and 5^2,m;wi,n, respectively. Here we use 
the vectorial notations Wi = {wi,i, ■ ■ ■ ,wi,p) and W2 — (^2,1, • ■ • ,W2,q). 

It is interesting to see what this involution means in terms of non-intersecting 
Brownian motions. As already discussed before, in the multiple Hermite case 
(jl.l8p - p.i9p the MOP are closely related to non-intersecting Brownian motions 
in the sense of Section 11.11 Then from the non-intersecting Brownian motions 
point of view, the above involution reduces to (cf. (|1.18p - (|1.19p ) 

t^l-t, P ^ q, {ai,ni) 'r^ {bj.rrij). 

This means that we interchange the role of starting and ending points and 
reverse the direction of time. In the te-plane this corresponds to a reflection of 
the Brownian motions with respect to the vertical line t — 1/2. 
In what follows we will denote 




(5.42) 
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Lemma 5.12. (See also fl [T^.J We have 

^^2,m;wi,n — "^^^1 ,n;W2,m"^ (5.43) 

where the superscript denotes the transposed inverse. 

The proof of Lemma 15.121 can be performed in a straightforward way. It 
suffices to check that both sides of (|5.43p satisfy the same RH problem, and 
next invoke the fact that the solution to this RH problem is unique; we omit 
the details. 

As an immediate corollary of Lemma 15.121 we can check what the involution 
does with the matrix Yi in (|5.ip . The Yi -matrices corresponding to the original 
and reversed Brownian motions turn out to be related as follows. 

Corollary 5.13. We have 

(^l)w2,m;wi,n = ^ ^(^1 ) wi ,n;W2 ■ (5.44) 

Put differently, if 

(Yi) wi,n;W2,m ^ 

denotes a partitioning with diagonal blocks of size p x p and q x q, respectively, 
then ^ ^ 

(Yi) W2,m;wi,n ( (J^^ — C"^ ' 

(5.46) 

Using Corollary 15.131 the compatibility relations in Section 15.3.21 may be 
immediately deduced from those in Section [5. 3. II 

5.4 Scalar product relations 

Now we discuss the scalar product relations induced by the compatibility con- 
ditions in Section [5.31 These relations are given by the following proposition. 

Proposition 5.14. (Scalar product relations) Let n and m with |n| = |m| he 
fixed and consider the partitioning (|5.24p for the matrix Yi = (Yi)n,m- Then 
we have the relations 

(C12C21),,, =t(l-t)^, (5.47) 

for all k {1, . . . ,p} , 

{C2iCu)^=t{l-t)'^, (5.48) 

for all I e {1, . . . ,q}, 

(Ci2C2i),,fe = -(1 - t){ak - a~,) (Cn),.,fc , (5.49) 
for all k,k ^ {1, . . . ,p} with k ^ k, and 

{C2iCi2)u = -t{bi - bj) {C22)u , (5-50) 
for all I, I £ {1, ...,(?} with I ^ I. 
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PROOF. If n = and m = then the solution Yq q of the RH problem (|1.20p - 
(I1.22P is upper triangular, so that C21 is the zero-matrix and the relations (|5.47p 
and (|5.48p hold in that case. For arbitrary n and m, these relations then follow 
from an induction argument based on ()5.3ip . (|5.38p . ()5.39p . and (15.4ip . 

The two relations ([QO]) and (1530]) are simply ([OS]) and (jOOl) . □ 

We call (I5.47p - (|5.50p scalar product relations, since these relations can be 
interpreted as giving the scalar products between row and column vectors of C12 
and C21. The relations ()5.47p and (|5.48p then involve rows and columns with 
the same index, while the relations (|5.49p and (|5.50p involve rows and columns 
with different indices. Moreover, note that these last two relations determine 
the off-diagonal entries of the diagonal blocks Cn and C22 in (|5.24p in terms of 
scalar products of the entries in the off-diagonal blocks C12 and C21 in (|5.24p . 

Let us illustrate the above scalar product relations for the 4x4 case (i.e., 
p — q ^ 2). Then the partitioning (j5.24p is, when written out in full, 

Cl,4\ 

C2,4 
C3,4 
C4,4/ 



(5.51) 

(5.52) 
(5.53) 
(5.54) 

In other words, the number of Brownian particles rii, n2 starting from the first 
and second starting point oi and 02, respectively, is (up to a common factor 
t{l — t)/N) directly expressed by the scalar products l|5.5ip . (|5.52p . while the 
number of Brownian particles mi, m2 arriving at the first and second endpoints 
hi and &2, respectively, is expressed by the scalar products (j5.53p . (|5.54p (up to 
the same common factor). 

The scalar product relations (|5.5ip - (|5.54p can be nicely expressed in terms 
of the Hadamard product matrix H of (|5.13p . Indeed, they express that the top 
rightmost 2x2 block of this matrix, i.e., 

\C2,3C3,2 C2,4C4,2/ 

has fixed row sums equal to t(l — i)^, t(l — i)^, and fixed column sums equal 
to t{l - t)^, t{l - t)^. It follows from this that the matrix ([535]) is fully 
determined by only one of its entries, since the other entries then follow from 
these row and column sum relations. This implies in particular the first three 
relations in Proposition 1 1.81 

In the case of general p and q, we have the following result. 



Cii C12 
C21 C22 





Cl,2 


Cl,3 


C2,l 


C2,2 


C2,3 


C3,l 


C3,2 


C3,3 


\C4,1 


C4,2 


C4,3 



Then (|5.47p and (|5.48p give us the four relations 



(C2,3 C2,4j 
(C3,l 03,2) 
(C4,l C4^2) 



C3,2 
C4,2 

Cl,3 
C2,3 

Cl,4 
C2,4 



t(l-i) 
t{l-t) 
t{l-t) 



ni 

TV' 

7^2 

iV' 

m.2 
TV 
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Proposition 5.15. (Row and column sum relations) Let n and m with |n| — 
|m| be fixed and consider the Hadamard product matrix H in (j5.12p . Then the 
top rightmost p x q block of the matrix H , i.e., the submatrix 



\^p,p-\-l^p-\-l,p • ■ • 



(5.56) 



has fixed row sums equal to t{l~t) 
to til-t)^,l^l,...,q. 



Cp,p+qCp+q,p/ 

. . ,p, and fixed column sums equal 



PROOF. These relations are equivalent with the scalar product relations 
(|5:i7l) and (15:151) . □ 



5.5 Spectral curve interpretation 



As an aside, let us discuss an equivalent formulation of Proposition l5.15l in terms 
of the so-called spectral curve. Let 



i^n.mCC, z) := det + -y„,m(z 



1. 



(5.57) 



where n = |n| = |m|. The zero set of the polynomial Pn,m(C, z) in (|5.57p defines 
an algebraic curve (the spectral curve). 

For example, when p = q = 2 we can use (|5.22|) - (j5.24p and (|5.57p to see that 



Pn,m(e,^)=detU/4 



N 



nt{l - t) 



det 



N (zh - Da ~Ci2 

nt{l -t)\ C21 Db 

z + (1 — t)ai ci,3 ci,4 \ 

^ - Z + (1 - t)a2 C2,3 C2,4 

-C3,l -C3,2 ^-tbi 

-C4,l -C4.2 |-t&2/ 

(5 



where C:= t{l-t)fC 

The spectral curve Pn,in{^, z) = in (|5.57p defines an algebraic curve which 
is of degree p + q in the variable ^. After resolution of singularities, this curve 
defines a Riemann surface which can be used as a tool in the steepest descent 
analysis of the RH problem. This is in fact a (p + (7)-sheeted Riemann surface 
where the ith sheet corresponds to the ith solution function ^(z) — ^i(z), i — 
l,...,p + q. 

Now it turns out [3j[7j[8j[T5l[3T] that an important role in the steepest descent 
analysis is played by the asymptotic expansions for z — > cx) of the functions ^i(z). 
It was observed for some special cases that the following asymptotic expansions 
hold. 

Proposition 5.16. (Asymptotic expansions of the spectral curve) Let p and q 
be arbitrary. Then when appropriately labeled, the p + q branches ^(z) = Ciiz), 
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i — 1, . . . ,p + q of the algebraic curve Pn,m{S,, z) = in (j5.57p behave as follows: 

nt(l — t) nt n z \z'^ J ^ 
, , , Nbi mil ^ / 1 \ 

as z — > CX3, /or any k S {1, . . . ,p} and I € {1, . . . , g}. 

Note in particular that the 1/z terms in the z —^ oo expansion, up to a 
factor ±1, are precisely the fractions of Brownian particles starting from the 
A:th starting point or arriving in the Ith ending point 6;, respectively. 

PROOF. We give the proof for the case p = q — 2; the proof for general p 
and q is similar. From (|5.58p . we see that the only way to have Pn,m{£., z) = 
when z — > cx) is that one of the diagonal entries in (|5.58|) vanishes as 2 — s- oo. 
This yields the required linear and constant terms in (|5.59p and (|5.60p . 

To find the 1/z terms in (|5.59p and (|5.60p . one can propose a term of the 
form c/z for some yet unknown c € C, and subsequently require the coefficient 
of the z term in the Laurent expansion of ()5.58p to vanish. For k = 1,2, we 
then getfor|fc=i(l-i)fa 

$fe(z) = z - (1 - t)ak ■ '— j 

which by ((?3T1) and ((53^ leads to ((53^)) . The 1/z term in ((530)) follows in a 
similar way from (|5.53p and (j5.54p . □ 

We end this subsection with a brief discussion. We described in Propo- 
sition 15.141 the scalar product relations induced by the Lax pair for multiple 
Hermite polynomials. In case where p = I these relations fully describe the 
off-diagonal entries of (FL)„^rn- Indeed, the block entries C12 and C21 in (|5.24p 
then have only one row and one column, respectively. The row and column sum 
relations (|5.47p and (|5.48p then give explicit expressions for the products CijCj^i, 
i < j. With a little bit more work one can also obtain explicit expressions for 
the individual entries Cij; this leads to the expressions in Similar remarks 
hold when q = I. 

In the case of general p and q, however, we are not able to find explicit 
expressions for the entries of (Yi)n,m anymore. In this case, the best that we 
could obtain is showing that the compatibility relations in Section 15.31 can be 
used to determine the off-diagonal entries of (yi)n+efc,m+e, in terms of the off- 
diagonal entries of (yi)n,m- This yields a recursive scheme for computing the 
off-diagonal entries of (FL)„^rn by induction on |n| — |m|. Unfortunately, we 
were not able to derive any useful asymptotic information from this recursive 
scheme; therefore we will not state these relations here. 

5.6 Completion of the proofs of Propositions 11.71 and 11.81 
5.6.1 Proof of Proposition [T77l 

As already noted above. Proposition 1 1 . 71 follows as a special case of Proposition 
\5.5\ except for the explicit form of the diagonal recurrence coefficients in the 
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(5.59) 
(5.60) 



first term in tlic riglit hand sides of each of (|1.39p - (|1.42p . The latter expressions 
are only valid in the multiple Hermite case. 

To derive these expressions for the diagonal recurrence coefficients, we will 
again use the compatibility conditions in Section 15.31 We proceed as follows. 
By evaluating the kth row of (|5.34|) . one obtains 

((1 - t)ak + tbiel + elCiiekel - e^Cii)Ci2e, = 0, 

which is 

(^{1 - t)ak + tbi + (Cii) ^ J {Ci2)kj - {CiiCi2)kj = 0, 
or equivalently 

From this expression, we obtain the relation 

in \ (n \ (^ +\ ^I^fe^^fc (Cii)fe,fe (Ci2)fe,i 

(^ii)fcfc~(^ii ^ [l - t)ak + tbi - 

' \ / k.k 



k,k (Cl2)fc,; 

Rewriting this in terms of the entries Cij, Cij of (yi)n,m and (yi)n+efc,m+e, , we 
have the equivalent expression 

^^fc^fc k^k r)4-l 

Ck,k - Ck,k ^ -t)ak +tbi ■ ■ , (5.61) 

Ck,p+l 

for the diagonal recurrence coefficients Ck,k ~ Ck.k- Note that the expressions in 
the first term in the right hand sides of each of p.39p - ()1.42p follows as a special 
case of (|5.6ip . This ends the proof of Proposition 11.71 

By the way, note that the right-hand side of (j5.6ip contains entries of (yi)n.m 
only, and so it is more convenient than (|5.14p which also contains an entry of 
(^2)n.m- In addition, (|5.6ip has only about half as many terms as (|5.14l) . 

Recall also that the left-hand side of ()5.6ip depends on I since Ck,k is an 
entry of the matrix (li)n+efc,m+ei • D 



5.6.2 Proof of Proposition [TTSl 

As already noted above. Proposition 11.81 follows from the more general row and 
column sum relations in Proposition 15.151 The only thing that remains is to 
prove (|L48l) . 

Assume again that p = q = 2. We will compute the determinant 

det C12 det C21 = det f ''^^^ ""^A det f ""^^^ """^'A (5.62) 

\C2,3 C2,4y C4,2y 
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first via 

det Ci2 det C21 = det(Ci2C2i) 

-dpf/^ ^(l-i)t -(l-t)(ai-a2)ci,2 

= - *)'^ + (1 - tfiai - a2)'ci,2C2.i, (5.63) 

where we made use of tlie scalar product relations (|5.47p and (|5.49p . 
We similarly compute 

det C12 det C21 = det(C2iCi2) 

= det M "i-^ 

^C44 €4,2/ \C2,3 C2,4 

i(6i-62)c3,4 
't{h-b2)ci,3 t(l-i)^ 



det , >, 7 f 

= - t)'^ + - fo2)'c3,4C4,3, (5.64) 

where we now used (|5.48p and (|5.50p . By equating (|5.63p and (|5.64p we obtain 
a relation between the off-diagonal recurrence coefficients ci,2C2,i and 03,404,3 

(1 - tfiai - a2)2ci,2C2a - t^bi - 52)'c3,4C4,3 = t'(l - t)'^^"~ 



iV2 

which for the special case ni = nii and 712 = 1112 reduces to (11.48^ . This ends 
the proof of Proposition [TTHl □ 
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